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Chapter  1 


INTRODUCTION 

The  electromagnetic  coupling  between  a  conducting  body  and  an 
aperture  is  an  important  problem.  One  of  the  studies  is  to 
investigate  the  response  of  protected  objects  behind  an 
aperture-perforated  metallic  screen.  The  power  transmitted  through 
the  aperture  from  a  source  on  a  nearby  object  which  is  behind  the 
screen  is  also  of  interest.  Situations  of  this  nature  arise  in  a 
number  of  applications,  e.g. ,  microwave  leakage  through  cracks  in 
shield  walls  of  electronic  equipment  and  electromagnetic  penetration 
into  vehicles,  aircraft,  and  ships.  This  penetration  is  through 
windows  or  seams  of  doors  to  cables  or  other  metallic  objects. 

In  the  past,  the  problem  of  electromagnetic  coupling  to  an 
infinitely  long  wire  through  an  electrically  small  aperture  [i]  -  [3l 
or  through  a  narrow  slot  [4]  has  been  studied.  The  problem  of  a  narrow 

slot  passing  by  a  finite-length  wire  without  loads  [5]  or  with  loads  f 6l 
has  also  been  investigated. 

In  this  work,  we  consider  more  general  problems  for  which  an 
aperture  in  an  infinite  conducting  plane  is  hacked  by  a  conducting 
body.  The  aperture  and  the  conducting  body  are  of  arbitrary  size  and 
shape.  A  moment  method  solution  is  developed  for  the  equivalent 
magnetic  current  in  the  aperture  and  the  electric  current  on  tbs 


-conducting  body.  The  solution  is  then  applied  to  problems  for  which 
the  conducting  body  is  a  wire  of  finite  length  (with  or  without  loads) 
or  of  infinite  length.  The  current  distributions  in  the  aperture  and 
on  the  wire  can  be  used  to  calculate  the  power  transmitted  through  the 
aperture  and  other  responses  of  the  wire. 

f 

I 

In  Chapter  2,  the  general  problem  of  a  conducting  body  behind  an 
|  aperture  is  considered.  The  equivalence  principle  [7,  Sec.  3-5], 

boundary  conditions,  and  moment  method  [s]  are  utilized  to  obtain  a 
pair  of  matrix  equations.  These  equations  are  used  to  solve  for  the 
current  distributions  in  the  aperture  and  on  the  conducting  body. 

In  Chapter  3,  we  specialize  the  problem  to  an  unloaded  wire  of 
finite  length  behind  an  aperture  of  arbitrary  size  and  shape. 
Triangular  patching  is  used  to  model  the  aperture.  Pulse  functions 
are  placed  on  the  wire.  Matrices  in  the  pair  of  matrix  equations 
developed  in  Chapter  2  are  evaluated.  Numerical  results  are  presented 
and  compared  with  those  in  C 5 "] . 

In  Chapter  4,  the  wire  considered  in  Chapter  3  is  extended  to  an 
infinitely  long  wire.  The  expansion  functions  and  matrices  developed 
in  Chapter  3  are  modified  to  include  the  effect  of  the  infinite  length 
of  the  wire.  Numerical  results  agree  very  well  with  the  available 
data  in  [l ]  -  [4]. 

In  Chapter  5,  the  problem  is  generalized  to  that  of  an 
arbitrarily  loaded  wire  behind  an  aperture.  We  modify  the  moment 
method  solution  developed  in  Chapter  4  to  include  the  reflection 
effect  from  the  loads  of  the  wire.  We  then  evaluate  the  power 


l 


transmitted  through  the  aperture.  The  solution  for  a  wire  with 
matched  loads  is  used  to  obtain  an  equivalent  circuit  of  the  aperture 
for  the  transmission  line  mode  on  an  arbitrarily  loaded  wire. 
Numerical  results  are  presented  and  compared  with  those  in  [4I  and 
[«]. 

He  summarize  and  conclude  this  study  in  Chapter  6. 
Recommendations  for  further  work  are  given.  Appendix  A  presents 
evaluations  of  some  surface  integrals  used  in  Chapter  3*  Fields 
due  to  outward  traveling  TEM  (transverse  electromagnetic)  currents  are 
derived  in  Appendix  H.  In  Appendix  C,  we  derive  some  integral 
formulations  used  in  Chapter  4.  An  equivalent  circuit  of  an  aperture 
for  the  transmission  line  mode  on  a  loaded  wire  is  derived  in  Appendix 
D. 

Figures  and  tables  are  at  the  end  of  each  chapter  and  appendix. 


Chapter  2 


FORMULATION  OF  THE  CORRAL  PROBLEM 


2.1 .  Problem  Specification 

The  general  configuration  is  given  in  Fig.  2.1.  An 

aperture-perforated  infinite  conducting  plane  of  zero  thickness 

separates  regions  a  and  b.  The  size  and  shape  of  the  aperture  are 

arbitrary.  In  region  a  (y  <  0),  there  are  impressed  sources 

(£iS,Mia).  In  region  b  (y  >  0) ,  there  are  an  arbitrarily  shaped  and 

il) 

sized  conducting  body  and  impressed  sources  (J1  f‘l  ).  The  material 

in  regions  a  and  b  is  loss-free  and  homogeneous  with  permeabilities 

(u  »u  )  and  permittivities  (e  ,e  ).  The  aperture  and  the  conducting 
a  b  a  b 

body  are  labeled  A  and  B,  respectively. 

In  this  chapter,  a  moment  method  solution  is  developed  to  obtain 


the  equivalent  magnetic  current  in  the  aperture  and  the  electric 
current  on  the  conducting  body. 
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2.2.  Formula  Derivation 

We  use  the  equivalence  principle  (7,  Sec.  3-5 1  to  divide  the 
problem  into  two  parts,  as  shown  in  Fig.  2.2.  The  fields  in  region  a 
remain  unchanged  if  the  aperture  is  closed  by  a  conducting  surface  and 
an  equivalent  magnetic  surface  current  sheet  M  is  placed  over  the 
aperture  region  where 


*  n.  x  over  A  (2-1) 

ri  is  the  unit  vector  normal  to  the  conducting  plane  and  points  toward 

region  b.  Ea  is  the  electric  field  over  the  aperture  in  region  a. 
“A 

This  equivalence  is  shown  in  Fig.  2.2(a).  Similarly,  the  fields  in 
region  b  remain  unchanged  if  the  aperture  is  closed  by  a  conducting 
surface  and  a  magnetic  current  sheet  j,s  placed  over  the  aperture 
region  where 


*  -  H  x  over  A  (2-2) 

E^  is  the  electric  field  over  A  in  region  b.  In  addition,  an  electric 
“"A 

current  J  is  induced  on  the  conducting  body  by  M  t  jib  f  and 
This  equivalence  is  shown  in  Fig.  2.2(b).  Therefore,  the  electric  and 
magnetic  fields  in  region  a  are 

(2-3) 

(2-4) 


£  -  E  (E  )  ♦  E 
Ha  -  f(!f)  ♦  Hia 


6 


Htre,  (Ea(Ma)t  lia(Ma))  is  the  field  from  ,  and  (E*S»  H**)  is  the 
short-circuit  field  (field  existing  in  the  presence  of  the  complete 
conducting  plane)  from  (£ia,  Mia).  Similarly,  the  fields  in  region  b 
are 

Rb  -  Eb(Mb)  +  Eb(j)  +  Eib  (2-5) 

Hb  «  Hh(Mb)  +  Hb( J)  +  Hib  (°-6) 

Here,  (^b(?lb),  H.b(Mb))  is  the  field  from  Hb,  (E<b(  J) ,  ]{b(j))  is  the 

field  from  and  (El  t  H^'3)  is  the  short-circuit  field  from  (J*b, 

M*b).  Note  that  all  fields  are  obtained  with  the  aperture  shorted. 

We  next  appply  boundary  conditions  to  this  problem.  They  are: 

(1)  the  tangential  electric  field  is  continuous  across  the  aperture  A, 

(2)  the  tangential  magnetic  field  is  continuous  across  A,  and  (3)  the 

tangential  electric  field  vanishes  on  the  surface  of  the  conducting 

body  Ba  The  first  condition  implies  that  n  x  E,  *  n  x  E.  in  (2-1)  and 

—  —A  —  —A 

a  b  a 

(2-2).  Therefore,  *  -M  .  From  now  on  we  use  M  to  represent  M_  and 

b 

-M  .  The  other  two  conditions  imply  that 

H®  -  hJ  on  A  (2-7) 

- t  - t 

_Eb  •  £  on  B  (2-8) 

where  t  denotes  the  component  tangent  to  A  or  B.  Substitution  of 
(2-4)  -  (2-6)  into  (2-7)  and  (2-8)  gives 
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■Hf(M)  -  JT(M)  ♦  ]r(j) 
— t  —  “t  “t  ~ 


,tia  ,.ib  . 
H  -  H  on  A 


Am)  -  e£(J)  =  j£b  on  B 
“t  - 1  — t 


(2-9) 


(2-10) 


Here,  the  linearity  of  the  operators  is  used  to  replace  H  (-M)  by 

— t  — 

-Hb(K)  and  jA-M)  by  -E^(M).  Equations  (2-9)  and  (2-10)  are  the 
t  L  u 

equations  to  solve  for  and  _J. 

We  next  use  the  moment  method  [s]  to  reduce  (2-9)  and  (2-10)  to 

matrix  equations.  For  this,  we  define  expansion  functions  {m  ; 

— n 

n*1 ,2, . . . ,NA}  over  A  and  jj  ;  n=1 , 2, . . . ,NB}  over  B  to  anproximate  M 

— n  — 

and  J_,'  respectively.  That  is, 


NA 

H  -  T  V  M 
“  L  n  “ n 
n-1 


(2-11) 


NB 

j  ■  y  i  j 

—  L  n  — r 
n«1 


(2-12) 


We  define  testing  functions  {Pi  ;  m*!  ,2, . . .  ,NA|  over  A  and  {5  ; 


—m 


m1 


•1,2,..., MB)  over  B  and  a  symmetric  product 


%  •  V 


» ■  n  *• 


.  F  dS 
-2 


(2-15) 


where  S  denotes  surface  A  or  B.  Now,  substituting  ( 2—1 1 )  and  (2-12) 

into  (2-9),  taking  the  symmetric  product  over  A  with  each  M  ,  and  then 

— m 

using  the  linearity  of  the  symmetric  product,  we  obtain  the  first 
matrix  equation.  Similarly,  substituting  (2-11)  and  (2-12)  into 
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(2-10),  taking  the  symnetric  product  over  B  with  each  j  ,  and  then 

m 

using  the  linearity  of  the  symmetric  product,  we  obtain  the  second 


matrix  equation.  The  two  matrix  equations  are 

[y*  .  y”1  v  .  Ft]  i  -  lla  -  lib 

(2-14) 

Ft]  v  *  [z]  I  -  vlb 

(2-15) 

Here, 

t**]  - 1  <-4  • 

(2-16) 

f'’”]  -  r  <-i  • 

(2-17) 

[t]  -  [  <?!  ,  «h(j  )>  ] 

— n  — t  ~ n  A  hAxHB 

(2-18) 

[t]  -  r  <j  ,  )>  l 

— n  — t  — n  B  HBxMA 

(2-19) 

[z]  -  r  <-j  ,  e^(j  )>j 

— m  — t  — n  3  NBxNB 

(2-20) 

■  r  <i  .  idVi  , 

-  — m  ~ t  A  NAxt 

(2-20 

"t-ib  r  /v  „ib.  -i 

1  '  £  <\  •  >»]-.*xl 

(2-22) 

^  m  f  <J  1 

L  h  3  HBxl 

(2-23) 

V  -  fv  ]  , 

'  n  NAxI 

(2-24) 

I 


Thus  (2-14)  and  (2-15)  are  the  matrix  equations  to  solve  for  the 

coefficients  (V  ,1  )  in  (2-1 1)  and  (2-12).  We  call  Ty3]  and  fYb]  the 
n  n 

admittance  matrices  for  regions  a  and  h,  respectively.  [t]  and  ^t] 
are  the  coupling  matrices,  [zl  the  impedance  matrix  ,  and  I*a  t  Tib, 

and  Vib  the  source  vectors.  V  and  I  are  the  unknown  coefficient 
vectors  to  be  determined. 


If  the  sources  in  region  a  were  to  radiate  in  the  homogeneous 

medium  characterized  by  (y  »£  )  and  if  the  sources  in  region  b  were 

a  a 

to  radiate  in  the  homogeneous  medium  characterized  by  (y.,e  ),  the 

°  0 

fields  in  (2-16)  -  (2-20)  would  be  given  by 


H1(jt_)  -  l«a,b 

—  — n  — n  n 


( 2-26) 


Eb(j  )  -  -juAb  -  V<J>b 


(2-27) 


£b(4)  -i- £ 

D 


(2-2R) 


Rb(M  )  -~VX  F* 
0 


(2-29) 


where 


F1  -  e  [[  M  (r ' )  G(k  r,  r1)  dS*  l-a,b 

— i n  l  j  j  “Ti  —  1 


(2-30) 


AD  -  u  J  ( r  * )  0(k  ,  r.  i')  dS' 


(2-31) 


10 


•i-JJ 


V-H  (rf) 

S,  c(1v  11  ^,)dS‘ 


l*a,b 


(2-32) 


V  f  ,  y’ . J  (l * ) 

4>„  ■  -  - ~n  —  G(k,  r,  r*)  dS1 

n  jj  Jueb  b*  -*  - 


(2-33) 


Vi 


0<kl’  i'  i'>  "  4tt | r -i '  t 


kl  ’  “Jel  *■! 


l*a,b 


l*a,b 


(2-34) 


(2-35) 


Here,  u»  is  the  angular  frequency.  All  fields  are  evaluated  in  a  half 
space  with  the  aperture  shorted. 

.  A  A 

If  a  Galerkin  solution  (i.e.,  M  -  M  ,  J  -J  )  is  used,  we  have 


[r1]  •  [r1]  i-«,b 


(2-36) 


[z]  -  [z] 


(2-37) 


[fl  -  -TtI 


(2-38) 


.SI 


where  ~  denotes  the  transpose  of  a  matrix.  This  will  save  us  some 
computations.  If  regions  a  and  b  are  filled  with  the  same  medium 
characterized  by  (y,  e),  we  can  remove  the  conducting  plane  and  use 
image  theory  [7,  Sec.  3-4]  to  evaluate  the  matrices  in  (2-16)  - 
(2-23).  They  become 


fra]  ■  frbl  -  2  F  <-L  .  ",(£)>,  L. 

— n  — t  — n  A  NAxNA 


(2-39) 


[T]  '  [  %  ■  V4  •  1 


NAxNB 


(2-40) 


(2-41) 


[T] 

[Z] 

lia 

I*b 


2  [  <3  ,  E  (H  )>  ] 

-m  -t-n  3  'NBxTJA 


f  <-L  ■  2.0,  •  2'  »*  1 


n  B  NBxNB 


2  f  <2»  •  4“>»  1 


A  NAxI 


2  r  <ft  »  Hjoh>.  L.  . 

Tn  — t  A  NAxi 


f  <J  »  E*ob+E*oh,>  ] 

1  ■ “m  -t  -t  B  JRBx1 


(2-42) 

(2- 43) 

(2-44) 

(2-45) 


Here,  H  (j  J*  )  and  E  (J  ,  J'  )  are  the  tangential  fielda  due  to 

— t  — n  —  n  — t  -n  ~~  n 

J  plus  its  image  J'  located  in  region  a.  H  (M  )  and  E  (M  )  are  the 

— n  —  n  — t-n— t-n 

i  o& 

tangential  fields  due  to  M  ,  H  is  the  tangential  field  due  to 

— n  — T- 

(J*a,  M^a).  H^°b  and  E*ob  are  the  tangential  fields  due  to  (j*b, 

M*b).  E^°b  is  the  tangential  field  due  to  the  image  of  (J*b,  n»*b). 

All  the  fields  are  evaluated  in  the  space  with  the  conducting  plane 

removed.  They  can  he  calculated  by  using  (2-26)  -  (2-35)  with 

superscripts  and  subscrints  (a,b)  dropped,  and  J  ieolaced  by  (J  ,.T'  ). 

“n  •  -n  -  n 


ka 


& 


a 


X  ~ 

M«»  ^0 

REGION  a 


CONDUCTING  PLANE 
(APERTURE  SHORTED) 


1 


3 


CONDUCTING  PLANE 


A 


Mb*  €l 


REGION  b 


Fi^.  2.2.  Equivalence  for  the  general  problem  .  (a)  Equivalence  for 

region  a.  (b)  Equivalence  for  region  b. 


sgwi 


V 

2 
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Chapter  3 


APPLICATION  TO  THE  PROBLRM  OF  AN  UNLOADED  WIRE  OF  FINITE 
LENGTH  BEHIND  AN  APERTURE  OF  ARBITRARY  SIZE  AND  SHAPE 


3.1  *  Problem  Specification 


In  this  chapter,  the  general  moment  method  solution  developed  in 

Chapter  2  is  specialized  for  the  problem  of  an  unloaded  wire  of  finite 

length  behind  an  aperture  of  arbitrary  size  and  shape.  The  problem 

configuration  is  shown  in  Fig.  3.1.  The  space  is  filled  with 

loss-free  homogeneous  medium  of  permeability  y  and  permittivity  e  .  A 

straight  thin  wire  of  radius  r  and  length  L  is  located  at  y*d  and 

B 

points  toward  the  z-direction.  A  plane  wave  is  incident  from  region  a 

at  an  angle  (e  ,d>  ).  For  this  excitation,  we  set  the  tangential 
o  o 

impressed  field  at  a  point  (x,0,z)  in  the  aperture  as 


(Hl08  u  ♦  I!1”8  u  1 

X  — x  *  — z 


(3-D 


where 


.  jk(xsin0  cos<j)  ♦zcosS  ) 

Hioa  -  (H^oacos6  cos<j>  -H*oasin<|>  )e  o  o  o 

X  0  O  0  <t>  0 


,ioa. 


Hioa  .  -Hioasine  e jk(x»in0rteos^Bco80o) 


(3-2) 

(3-3) 
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Here,  H*oa  and  H*°a  are  constants,  and  u  and  u  are  unit  vectors  in 

e  <t>  -x  -z 

the  x  and  the  z  directions  of  the  rectangular  coordinate  system.  Ho 
impressed  sources  exist  in  region  b. 

An  appropriate  set  of  expansion  functions  and  testing  procedures 
are  needed  to  solve  for  the  currents  in  the  aperture  and  on  the  wire. 
In  this  chapter  and  throughout  the  following  chapters,  a  Galerkin 
solution  is  utilized  and,  for  simplicity,  the  same  medium  is  assumed 
in  regions  a  and  b.  Therefore,  we  only  need  to  specify  the  expansion 
functions  and  use  (2-39)  -  (2-45)  to  evaluate  the  matrices. 


3.2.  Magnetic  Current  Expansion  Functions  in  the  Aperture 

In  this  section,  we  model  the  aperture  by  planar  triangular 
patching  which  has  been  used  to  model  a  surface  in  [9]  -  [13]«  The 
advantages  of  the  triangular  patching  are  :  (1)  the  ability  to 

conform  closely  to  an  arbitrarily  shaped  aperture,  (2)  the  flexibility 
of  having  greater  densities  on  those  portioas  of  the  aperture  where 
more  resolution  is  desired,  e.g.,  when  the  edge  effect  is  concerned, 
and  (3)  the  scheme  can  be  easily  implemented  on  a  computer.  The 
triangular  patch  scheme  for  the  aperture  is  explained  as  follows. 

First,  assume  that  a  suitable  triangulation  i3  found  to  closely 
approximate  the  aperture  region.  The  triangulation  is  defined  by  sets 
of  faces  (patches),  edges,  and  nodes,  such  as  shown  in  Fig.  3*2.  We 
number  the  nodes  and  edges  and  specify  the  orientations  of  edges  by 


arrows  (e.g. ,  edge  1  ia  from  node  3  to  node  4).  We  next  define  the 
orientation  of  each  face  to  be  its  normal  direction,  i.e. ,  in  the 
y-direction  of  Pig.  3.1.  As  shown  in  Fig.  3.3,  the  current  reference 
direction  across  an  internal  (non-boundary)  edge  n  is  defined  to  he 
the  direction  of  the  cross  product  of  the  edge  orientation  and  the 
face  orientation.  The  conjoined  triangles  associated  with  edge  n  are 

4  -  4 

denoted  by  T  and  T  ,  with  the  current  reference  direction  from  T  to 
n  n  n 


Next,  we  introduce  local  position  vectors  associated  with  edge  n, 

as  shown  in  Fig.  3.4.  Any  point  in  triangle  T+  can  be  designated  by  a 

n 

local  position  vector  defined  with  respect  to  the  free  node  (node 

not  on  edge  n)  of  T+.  Similarly,  any  point  in  T~  can  be  designated 

n  n 

by  £_n>  Fote  that  is  defined  away  from  the  free  node  of  Tfl,  and 

is  toward  the  free  node  of  T".  r  is  the  global  position  vector  of  a 

n  — 

point  in  the  aperture. 


Now,  we  define  the  expansion  functions  in  the  aperture,  which 
were  originally  proposed  in  f  1 4.1 .  For  each  internal  edge  n,  the 
magnetic  current  expansion  function  is  defined  a3 


r  in  T 


131 


>I(r) 

J1  — 


r  in  T' 


(3-4) 


elsewhere 
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Here,  1  is  the  length  of  edge  n,  and  A  and  A  are  the  areas  of 
n  n  n 

+  - 

conjoined  triangles  and  T  . 


In  the  following,  we  discuss  some  nroperties  of  M  defined  in 

— n 

(3-4),  and  ensure  that  they  are  uniquely  suited  to  be  the  magnetic 
current  expansion  functions  in  the  aperture. 


(t)  M  has  no  component  normal  to  the  boundary  edges  (excluding  the 
— n 

common  edge)  of  the  conjoined  triangles  and  T~.  Hence,  no 

n  n 

line  charges  exist  along  these  boundaiy  edges.  In  addition,  the 
component  of  M  normal  to  edge  n  is  the  constant  1  (see 

“T1 

Fig.  3.5)  and  is  continuous  across  edge  n.  This  implies  that 
edge  n  is  also  free  of  line  charge. 

(2)  The  surface  magnetic  charge  density  i3  defined  as 


7*M 

m  *  ■ 

..  n 

n 

-jU) 

+1 

_ _ n_ 

r  in  T“ 

—  n 

-  * 

0 

elsewhere 

(3-5) 

Therefor 

k 

e,  the 

net  charge  in  the  conjoined 

triangles  is  zero. 
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(3)  The  set  of  functions  M  defined  for  all  the  internal  edges  is 

— n 

sufficient  to  represent  the  magnetic  current  anywhere  in  the 
aperture.  Any  current  in  the  aperture  can  be  well  approximated 
by  a  superposition  of  these  expansion  functions.  At  any 
boundary  edge  of  the  aperture,  the  normal  component  of  magnetic 
current  must  vanish,  while  the  tangential  component  can  be 
represented  by  a  linear  combination  of  the  expansion  functions 
associated  with  the  two  internal  edges  of  the  pertinent 
triangle.  Therefore,  we  need  not  and  should  not  define 
expansion  functions  for  the  boundary  edges. 

Thus,  the  total  number  of  expansion  functions  over  A,  NA  in 

(2-11),  is  equal  to  the  total  number  of  internal  edges  of  the 

aperture.  In  addition,  since  the  normal  connonent  of  M  across  edge  n 

— n 

is  1 ,  V  in  (2-11)  is  interpreted  as  the  normal  component  of  the 
n 

magnetic  current  density  crossing  edge  n. 


3«3«  Electric  Cur  rent  Expansion  Functions  on  the  Wire 

We  make  the  following  approximations  for  the  current  on  the  wire  : 
(l)  The  current  is  assumed  to  flow  only  in  the  axial  direction  of 
the  wire,  (2)  the  current  is  approximated  by  a  filament  of  current  on 
the  axis  of  the  wire  and  depends  only  on  the  axial  length  variable, 
and  (3)  the  boundary  condition  (2-10)  is  applied  only  to  the  axial 
component  of  the  electric  field  at  the  surface  of  the  wire.  For 


simplicity,  the  wire  is  equally  divided  into  *JB  subsections  of  length 


2a  .  Pulse  functions  are  used  over  each  subsection,  "’he  expansion 
B 


functions  for  the  current  J  on  the  wire  are  then  defined  as 


J.  *  u.  P(z-z  )  n-1 ,2 

n  z  n 


,NB 


a  <  z  <  =■  - 
aB  ~  ~2 


°3 


P(z-z  ) 
n 


,  z  <  z  <  z 
n-  —  ~  n+ 


,  elsewhere 


(5-6) 


(3-7) 


where  z  ,  z  ,  and  z  are  the  z-eoordinates  of  the  starting  point, 
n-  n  n+ 

midpoint, and  termination  point  of  the  nth  subsection.  The  subsections 
are  enumerated  from  z  -  -L/2  *a^  to  L/2  -ag*  Note  that  by  (3-6)  the 

current  is  ensured  to  be  zero  at  the  ends  of  the  wire. 

j 

I 


3.4.  Matrix  Evaluation 


(A)  Evaluation  of  Admittance  Matrix  [Ya  _+  Y^l 

According  to  the  definition,  the  admittance  matrix  [Ya+yh]  is 
independent  of  excitations  applied  to  the  system  and  can  be  evaluated 
as  if  the  wire  were  not  present.  Therefore  the  matrix  is  the  same  as 
that  evaluated  in  1*131  •  We  summarize  the  evaluation  a3  follows. 
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Substituting  (2-26),  (2-30),  (2-32),  (3-4),  and  (3-5)  into 

(2-39),  we  have 


Ya  +Yb  -  4  jo)  f  [  M  *F  dS  +  4  f  K  .V4>  dS 

mnnn  J  J  — m  ~n  j  J  ~ m  n 


4  j  to  M  •  F  dS  -  4 


If  4>  V*  M  dS 

Jj  n  t 


1  S  «P  rt  1  p  *F 

4jo>r  ff  -n--Va-ds  ♦  m— w  -zJX.  dS  1 

JJ  2A  JJ  2A 


1  $ 

[ff -vas  -  r 

JJ  A_  JJ 


2 jw  1  TpC+*F  (rc+)  +  0C"*F  (tc')l 
n  '“in  — n  -to  -to  — n  — m 


41  0  (rc+)  -  $  (r®“)] 
m  n  -to  n  —to 


(3 -ft) 


where 


F  (*!?>  "  e  ff  «  (r ' )  C(k,  r®4,  r’)  dS* 
— n  —to  +J  J  _  — n  —  -to  — 


T  +T 
n  n 


(3-9) 


*  (iff)  “-n*  ff  mn(jL*)  G(k.  x?.  t')  dS* 
n-ro  y  +J I  n  to  — 

T  +TO- 
n  n 


(3-10) 


and  _r°+  are,  respectively,  the  local  and  global  position  vectors 

of  the  centroids  of  1 T.  The  integrals  are  approximated  by  sampling 

m 


the  integrands  at  the  centroids  in  (3-B). 


To  evaluate  P  (r0-)  and  $  (r0-),  we  proceed  face  by  face  for  each 
— n  -m  n  — ® 

triangle.  As  shown  in  Fig.  3.6*  T  is  an  observation  triangle  and  T 

p  q 

a  source  triangle.  Nodes  of  T  are  designated  by  1 ,  2,  and  3,  and 

q 

lengths  of  edges  of  T  are  1  ,  1  .  and  1  .  The  local  position  vectors 

q  12  3 

of  a  source  point  in  T  ,  p  ,  6  ,  and  3  *  divide  T  into  three  subareas 

q  -1  -2  -3  q 

A1 .  A^,  and  A^.  The  global  position  vectors  of  the  source  point  and 

the  three  nodes  of  T  are  r',  r  ,  r  ,  and  r  ,  resDectively.  We  now 

q  -  -1  -2  ~3 

define  area  coordinates  (£,n*£)  as 


(3-n) 


(3-12) 


-  5  -  n 


(3-13) 


where  A  is  the  area  of  T  .  The  relationships  anong  the  area 

q  q 

coordinates  and  position  vectors  are 


(Ll  m  ±(L'  -  O  i  -  1  ,  2,  3 


(3-14) 


L'  “  Ci,  +  Cl,  +  nr. 


(3-15) 


In  (3-14),  the  positive  sign  is  used  if  the  current  reference 


direction  of  edge  i  is  away  fron  T  ,  and  the  negative  sign  is  used  if 

q 

the  current  reference  direction  is  toward  T  .  Note,  for  simplicity, 

q 

+  + 

superscripts  ♦  of  p-  and  T  are  drooped. 

-i  q 
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Using  the  area  coordinates,  we  can  transform  any  surface  integral 
over  T  into  a  double  line  integral  by  the  following  formula, 

q 


jj  f(r_’)dS'  -  2A^  J  J  f((l-?-n)r_i^r2+nr^)d?(n 


r\  *0  K  B0 


(3-16) 


Therefore,  the  electric  vector  potential  and  magnetic  scalar 

potential  $  at  the  centroid  rc  of  triangle  T  due  to  edge  i  of 

l  — p  p 

source  triangle  T  can  be  evaluated  as  follows. 

q 


3y  (3-4),  (3-9),  and  (3-1-1)  -  (3-16),  we  have  the  vector 

potential 


FP9  -e  ff  K  (r ' )  G(k,  rc,  r’)  dS’ 

-a  J)  ~ i  ~  T  ~ 


-ilff 

2A„  JJ 


p.  G(k,  r  ,  r’)  dS* 
~1  ""P  ““ 


± 1  G  fr 

~~  f(l-C-n)r  +?r  +Tir  -r  ]G(k,q  )dS' 
2Aq  JJ  -1  -2  -3  “I  p 

T 

q 


±1i£f  >-  ,3  (3-17) 


where 


i  i  -n 


|  j  G(k,Rp)d?dn 


(3-18) 


n-o  C-o 


€  nio  £*o 


J  J  CG(lc,R  )dcdn 


(5-19) 


ipq  «  j  J  nG(k,Rp)d£;dn 

^  n-o  £>o 

-jkR 

V ' 


(3-20) 


(3-21) 


Rp  ■  ll!  -  0  -  K  -  n)jL,  -  S  L2  -  9  i3l 


(3-22) 


Equations  (3-18)  -  (3-20)  are  evaluated  in  Appendix  A.  Similarly,  by 
(3-5)  and  (3-10)  -  (3-i5),  we  have  the  scalar  potential 


*pq  fk(r’)  G(k,  r°,  r/)dS* 
^  ii  i  —  —v  ~ 


j-  |J  0(k.  i°.  i')  SS1 


f-1.  IM  1-1 ,2,5 
Juy  i 


(3-23) 


Again,  if  the  current  reference  direction  of  edge  i  is  away  from  T  , 

q 

the  positive  sign  in  (3-17)  and  the  negative  sign  in  (3-23)  are  used. 
If  the  direction  is  toward  T  ,  the  negative  sign  in  (3-17)  and  the 

q 

positive  sign  in  (3-23)  are  used.  Note  that  (3-17)  and  (3-23)  are 
defined  only  for  the  edge  which  is  not  a  boundary  edge  of  the 
aperture.  If  edge  i  is  a  boundary  edge,  F?q  -  <5>  ^q  *  0. 
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.  4 

t 


Finally,  F  (rc  )  and  $  (ic  )  in  (3-9)  and  (3-10)  can  he  obtained 
— n  — m  n  “m 

fron  fP**  and  in  (3-17)  and  (3-23)  by  transferring  the  local  index 

i  (i  *  1,  2,  3)  to  its  corresponding  global  index  n  (n  *  1,  2,..., 

NA).  Similarly,  the  quantity  lj^j1  iri  (3-9)  can  be  evaluated  by  using 

face  to  face  procedure  and  local  index.  The  evaluation  of  the 

admittance  matrix  is  now  complete. 


(B)  Evaluation  of  Impedance  Matrix  [ z] 


The  tangential  electric  field  at  an  observation  point  along  the 

surface  of  the  wire  due  to  the  axial  current  J  at  (0,d,z’)  dIus  its 

n 

image  J'  “-J  at  (0,-d,z')  can  be  obtained  by  (2-27),  (2-31),  (2-33), 

-  n  — n 

and  (2-34).  It  is 


E  (J  ,J'  )- 
z  — n  —  n 


47T 


♦  — 1—  j- 

4Tr)jo)  e  8Z 


R* 


)dz  * 


dJ  (z*)  -jkP 

_JJ -  (S - 

dz'  R 

B 


)  dz ' 


n*1  ,2, . . .  ,TIB 


(3-24) 


Here, 


z')2  ♦  (2d)2 


(3-29) 


R*  - 


(3-26) 
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v> 

& 


\y 

l  V 


aze 

v: 

to 

$ 

& 

and 

impedance  matrix  elements 


Z  *  jtoViAl  Al  T  tp( a  ,  z  ,i  )  — x{j( z  ,  z  ,2d)] 
mn  m  n  “  —  m  n  B  ~mn 


+  -nb[i(z  *'z  .  * ^t> ^  "  !t(z  **z  ,.»2d) 

joje  — *  m+  n+  B  m+  n+ 


\p(z  ,z  ,r  )  +  Ip(  z  ,z  ,  2d) 
—  m+  n-  B  —  m+  n- 


it>(z  ,z  ,r  )  +  »|/(z  ,z  ,2d) 
—  m-  n+  B  —  m-  n+ 


(3-27) 


+  ty(z  ,z  ,r  )  -  ^(z  t z  ,2d) ]  mf  n*1,2,...,NB 
—  m-  n-  B  n-  n- 


where  Al  -Al  *2a  and 
m  n  B 


Al 


2  ♦  —  p- .  ■  . 

n  2  -jkjz  -z’)2+P2 

r  m 


^(z  ,  z  ,  P)  -  [  f-  -- 

m  n  4ttA1  J  ,  .  \ 

n  Al  J^VZ  } 


2-p2 


dz’ 


(3-28) 


z - i 

n  2 


has  been  evaluated  in  [l5]» 
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A 

(C)  Evaluation  of  Coupling  Matrices  [t]  and  PtJ 


The  tangential  magnetic  field  at  a  point  (x,0,z)  in  the  aperture 
due  to  a  z-  directed  current  J  on  the  wire  at  y*d  plus  its  image 


H  (J  ,  J'  )  -  2H  (J  )  u 
~t  ~n  ”  n  x  ~ n  ~x 


f  "  w  H  (J  )  u 


(3-29) 


Substitution  of  (3-4)  and  (3-29)  into  (2-40)  gives 


.  rt  p+*u  HA(J  )  .  rr  p  ’u 

-di  ds  .i  ^ 

m  a  J  2  2  A  "  2 

b  4  Jx  +d  m  _  J  x  ♦ 


3"*u  h,(j  ) 

-m  —x  9  n 


t  ffw  *r>  , 

ft*  m+)2+d?  ftxpW 


(3-30) 


Here,  £c  and  x°  are  the  x-components  of  oosition  vectors  Pc+  and 
ra  m  ■■ 

c±  .  /  c±  \ 

,  respectively.  H^( J^;  r^  )  is  the  (^-component  of  the  magnetic 
field  due  to  J  at  the  centroid  rc_. 


The  magnetic  field  H  (j  ;  rC±)  due  to  J  defined  in  (3-6)  is 

<J>  — n  — m  — n 


H.(j  ; 

4>  — n  y3p 


(3-30 


A1 

n 


(3-32) 


is  the  distance  between  the  centroid  associated  with  edge  m  of  the 
aperture  and  the  midpoint  of  subsection  n  of  the  wire.  zc±  is  the  z- 

HI 

component  of  x?~ .  Substituting  (3-31)  into  (3-30),  we  obtain 


dl  A1 


™  nj  „c+r  jk  *  1  i 

^  "  1Z7  (.Djl 


♦  fj-  t_£L_.  '  1  e 


a 

— ,  1  e  ""] 


m"> • • « fNA  n-1 f • • • ^NB 


(3-33) 


Since  Galerkin's  solution  is  used,  [t]  can  be  obtained  (t]  »  -  [t1. 


! 


/  .  ->-ia  -*-ib  -►ib 

(Dj  Evaluation  of  Source  Vectors  I ,  I ,  and  V _ 


Since  there  is  no  excitation  from  region  b  of  Fig.  34 
(2-44)  and  (2-45)  give 


-*ib  -►ib 
I  -  V  -  0 


Substitution  of  (3-0  -  (3-4)  into  (2-43)  gives 


4a 


I  -  1 


m  m 


[-7  ||  C^oa  ds*-L:  II  £-4“‘sl 


m 


-  1  [x  H10a(rc  )  ♦  z°  Hloa(r c+) 
m  m  x  — m  m  7.  — m 


♦  «c"Hloa(rc~)  ♦  «c"Hioa(rc’)  ] 
m  x  — m  m  z  — m 


m-1 ,2, ... ,NA 


Here, 


H^0a(r^~)  «  (Ha°acos0  cos(J)  -  H^oasin<f>  ) 
x  D  009  0 


,1k(xc+sin0  cos«j)  +  zc~cos0  ) 


c± 


m  0  0  1:1  0 


Ana  /.♦  -tno  jk(xG±sin0  cos 41  +ZC±COS0  ) 

Hioa(i.ci)  -  -  HnOasin0  e  m  00m  o 
z  m  o  o 


.c± 


and  Zm  ar®  the  ®“comP°n®nt8  of  pc±  and  r0* , 


(3-34 


(3-351 


(3-36: 

(3-37; 


respectively. 


£ 


By  solving  (2-14)  and  (2-15)  with  matrices  evaluated  above,  we 


1 

1 


can  obtain  the  current  distributions  in  the  aperture  and  on  the  wire. 


3.5.  Numerical  Results 


In  this  section,  numerical  results  are  presented  for  the  magnetic 
current  distribution  in  the  aperture  and  the  electric  current 
distribution  on  the  wire.  Numerical  results  for  narrow  3lots  are 
available  in  the  literature  for  the  problem  considered  in  this 
chapter.  We  therefore  perform  computations  for  a  narrow  slot  backed 
by  an  unloaded  wire  of  finite  length  with  plane  wave  incidence.  As  we 
will  see,  our  results  agree  very  well  with  those  in  [5]. 


& 

& 


The  slot  has  length  L  and  width  W  and  is  centered  at  (x  ,0,0). 

a  a  c 

The  wire  has  length  L  and  radius  r  and  is  located  at  y“d.  A  plane 

B 

wave  is  normally  incident  with  H^oa“l  ampere/meter.  Figure  3.7  shows 
the  triangulation  used  for  the  slot,  where  NA,  TJNODE,  HEDGE,  and  NFACE 
are  the  numbers  of  unknowns,  nodes,  edges,  and  faces  of  the  slot. 


respect ively. 


1 


For  the  case  of  L  /X  “0.5,  W  /X  “0.05,  r  /X  “0.001,  and  L 
a  a  B 

/X-0.5,  Fig.  3.B  shows  the  comparison  of  our  results  for  the  total 
x-directed  magnetic  current  M  in  the  3lot  with  those  in  f 5 1 •  Figure 


I? 

£ 


3.9  shows  the  comparison  of  the  electric  current  I  on  the  wire.  They 
both  have  very  good  agreement.  Figure  3.10  shows  the  electric  current 


on  the  wire  for  the  case  of  L  /\  -0.5,  W  /X  “0.05,  x  «0.,  d /X  “0.1 25, 


3 

I  vs 


1 


•  *.•*•*•'  ■ .  ■  ,  *  .*•  • .  *  •'*  »**  • 


¥ 
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r  /X“0.001 ,  and  L/  X=i .0.  Figure  3.1 1 
B 

the  wire  for  the  case  of  L  /  X^'.O,  V  j 

a  * 

r  /  X“^»001 ,  and  L  /  A=1  *0.  A/»nin,  our  results  for  these  cases  aeree 
very  well  with  those  in  [5].  Here,  X=*1  meter  is  chosen  as  the 
wavelength  in  the  computation.  Hote,  the  opposite  siyn  of  electric 
current  on  the  wire  in  our  results  is  due  to  the  choice  of  a  different 
coordinate  system  from  that  in  f 5 1 • 


shows  the  electric  current  on 
X=0.05,  xj  X-0.25,  d/  X=0.25, 


A  * 

If. 


Pig.  3.1*  An  unloaded  wire  of  finite  length  behind  an  aperture. 
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@^~NODE 


ORIENTATION  OF  EDGE 


Fig.  3*2.  Trian/nilation  example. 


EDGE  ORIENTATION 


FACE  ORIENTATION 
(NORMAL  TO  THE  FACE 
POINTS  OUT  OF  THE 
PAPER) 


CURRENT  REFERENCE  DIRECTION 


Fig.  5.3.  Relationship  among  face  orientation,  edge  orientation 
current,  refeience  direction  across  edge  n. 


EDGE  ORIENTATION 


Fig.  3.4.  Local  position  vectors 


— n 


and  p-  associated  with  edge  n,  and 
— n 


global  position  vector  r. 


Fig.  3.5.  Normal  components  of  3  and  0  at  edge  n. 
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Fip.  3»7.  Triangulation  of  a  slot  centered  at  (x  ,0,0)  for  NA*19, 

c 


NNODK-22,  NRDGK-41,  and  NFACR-20. 


- OUR  RESULTS (NA*I9,  NB*9) 

XA  RESULTS  IN  [5] 

Fig.  3.8.  The  total  x-directed  magnetic  current  in  a  slot  backed  by  an 

unloaded  wire  of  finite  length  with  normal  plane  wave 

incidence.  L  /A*0.5»  W  /^*0.05»  x  dM*0.25» 

a  a  c 

r  /x-O.OOt,  L/X-0.5,  and  Hioa-1  ampere/meter.  (X-1  meter) 


-  OUR  RESULTS  (NA* 19,  NB«9) 

X  A  RESULTS  IN  [5] 

Fig.  3-9*  The  electric  current  on  an  unloaded  wire  of  finite  length 

passing  by  a  slot  with  normal  plane  wave  incidence. 

I  /X-0.5,  W  /A-0.05,  x  -0.,  d/X-0.25,  r „/X«0.001 ,  L/X-0.5, 

a  a  c  B 

and  H*oa*1  ampere/meter.  (X*1  meter) 


-  OUR  RESULTS  (NA«I9,  NB*I9) 

x  A  RESULTS  IN  [5] 


Fig.  5.10.  The  electric  current  on  an  unloaded  wire  of  finite  length 
passing  by  a  slot  with  normal  plane  wave  incidence. 
L  /X-0.5,  W  A-0.05,  x  -0.,  d/X-0.125,  1  A-0.001 , 

a  a  c  u 

and  Hioa-1  ampere/meter.  (X  “1  meter) 
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-  OUR  RESULTS  (NA«I9,  NB«24) 

xA  RESULTS  IN  [5] 


Pig.  3.11.  The  electric  current  on  an  unloaded  wire  of  finite  length 

passing  by  a  slot  with  normal  plane  wave  incidence. 

L  A-1.0,  W  /A-0.05,  x  /A-0.25,  d/A«0.25,  r  /A-0.001 , 

a  a  c  E 

L/A*1 .0,  and  H*oa-1  ampere/neter.  (A  *1  meter) 
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Chapter  4 

APPLICATION  TO  THE  PROBLEM  OF  AN  INFINITELY  LONG 
WIRE  BEHIND  AN  APERTURE  OF  ARBITRARY  SIZE  AND  SHAPE 


4.1 •  Problem  Specification 


In  this  chapter,  the  wire  considered  in  Chapter  3  is  extended  to 
an  infinitely  long  wire.  The  geometry  of  the  wire  and  an  aperture  of 
arbitrary  size  and  shape  is  shown  in  Fig.  4.1.  The  excitation  is  a 
plane  wave  incident  from  region  a  of  Fig.  4.1 ,  and  is  defined  in 
(3-1).  The  expansion  functions  and  the  matrices  developed  in  Chapter  3 
are  modified  to  include  the  effect  of  the  infinitely  long  wire.  For 
this,  two  exponential  wave  functions,  in  addition  to  the  pulse 
functions,  are  added  to  the  expansion  functions  on  the  wire.  The 
objective  is  to  obtain  the  current  distributions  in  the  aperture  and 
on  the  wire. 


4.2.  Expansion  Functions 

Since  the  aperture  is  arbitrarily  sized  and  shaped,  the  expansion 
functions  for  the  aperture  remain  the  same  as  those  developed  in 
Chapter  3*  It  is  noted  that  a  transmission  line  without  reflections 


is  formed  by  the  presence  of  the  infinitely  Ions  wire  and  the 
conducting  plane.  Therefore,  the  electric  current  on  the  wire 
consists  of  two  outward  traveling  TEM  currents  plus  evanescent 
currents  (higher  order  modes)  existing  in  a  finite  region  near  the 
aperture.  We  use  pulse  functions  to  represent  the  total  current 
(evanescent  currents  plus  outward  traveling  TEM  currents)  in  this 
finite  region.  We  use  two  exponential  wave  functions  to  represent  the 
outward  traveling  TEM  currents  outside  the  finite  region.  This  finite 
region  must  be  greater  than  the  dimension  of  the  aperture  in  the 
z-direction.  Therefore,  the  electric  current  expansion  functions  on 
the  wire  are 


u  e 


jkz 


oo  <  !  (  • 


(4-1) 


J  -  u  P(z-z  ) 
~~n  ~z  n 


|z|  1  j  ,  n-2,3 . NB-i 


(4-2) 


<J  *  u  e 
—MB  ~z 


-jkz 


z  <  00 


(4-3) 


Here,  P  is  the  pulse  function  defined  in  (3-7)  and  L  is  the  length  of 
the  finite  region.  Each  subsection  in  this  finite  region  is  of  length 
2ctg«  Note  that  in  order  to  include  the  outward  traveling  TEM 
currents,  the  indexes  in  (4-2)  are  slightly  changed  from  those  of 

(3-6). 
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4.3.  MATRIX  EVALUATION 


It  is  noted  that  the  magnetic  current  expansion  functions  for  the 

aperture  and  the  electric  current  expansion  functions  for  the  finite 

region  of  the  wire  are  the  same  as  those  in  Chapter  3.  Therefore,  the 

corresponding  elements  of  the  matrices  [y®+y^],  "lia,  an(j  are 

identical  to  those  in  (3-8),  (3-34),  and  (3-35).  The  elements  Z  for 

mn 

m,  n  •  2,  3 . NB-1  are  defined  in  (3-27),  while  the  elements  T  ■ 

mn 

-  T  for  m  ■  1,  2,...,  NA  and  n  ■  2,  3,  ...,  NB-1  are  defined  in 

nm 

(3-33). 


We  now  evaluate  the  matrix  elements  Z  associated  with  the 

mn 

outward  traveling  TEM  currents  on  the  wire.  Ry  ( B— 1 7)  and  (B-19),  the 
z-componenta  of  electric  fields  due  to  plus  its  image  £’ and 
4b  plus  It.  in*.  «re 

'•1k2  1|t  ..t  -jkr' 

rQ  (r0} 


.  [JL  -  Z*2  lc"3kr  < 

2  /_ i -\3 


(r’o>  ^ 


(4-4) 


f^ir  -4- 

4irjaje  ♦  (r*)2  ,  ^3 

r0  '0'  'r0/ 


t'o  (r’o)2  (r^)5 


(4-5) 


where 
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+ 


(4-6) 


r 


»  + 

0 


J<«*  j)2  *  (2d)2 


(4-7) 


However  it  is  noted  that  the  charge  densities  associated  with 
and  J  _  are 


-NB 

■jwdz 

-Mike3**  -  J  )e'Jk2l 

(4-9) 

dJ 

_ NIL* 

-jwdz 

(4-9) 

The  z-components  of  electric  fields  due  to  point  charges 

[+6(z?L/2)e"*)kL/2]/(-ja))  plus  their  images  are  the  terms  involving 

(Iq)  2,  (r*)  (r'*)  2|  and  (r’^)  ^  in  (4-4)  and  (4-6).  To  assure 

the  continuity  of  the  total  current  on  the  wire,  we  first  replace 

delta  functions  6(z+L/2)  in  (4-9)  and  (4-9)  by  pulse  functions 

P(z+L/2)/A1a  over  a  small  region  Al  .  Then  we  replace  the  E  due  to 

the  point  charges  by  the  E  due  to  these  pulse  functions.  Al  is 

z  0 

chosen,  for  simplicity,  to  be  equal  to  2a  .  Thus,  (4-4)  and  (4-6) 

B 

become 


-jkr~  -jkr'" 

E  (J.,J*  )  - (~ - - - )  ♦  E  (P-) 

z  — 1  —  1  4moe  -  i-  z 

r0  0 


-.1k§  -Jkrj  -Jk,** 

Ez^NB  -  l?B^  "  4ttw£  ^  ♦  T^~^  +  I,'z^P  ^ 

ro  1  0 


(4-10) 


(4-11) 
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Here, 


-frj 

E  (p*)  -  — - 

V  '  4irjtoeAl 


±e _ 1_ 

4TTja3eAL0  32  J 


a  f  t  _"JkR  -jkR’ 

ij  -* — 


R 


-jkr±  -jkr’± 

P(2')(^—  -  -§ -  W 

r±  r ' ± 


(4-12) 


to 


is  obtained  from  the  second  term  of  (3-24).  R  and  R'  are  defined  in 
(3-25)  and  (3-26),  and 


f. 

te 


2  ”  '  +IB 


■*  - 


(4-13) 


(4-14) 


By  (2-42),  (4-1),  (4-10),  (4-12),  (C-3),  and  (C-4),  we  have  . 


-jkL 

Zn  -  *  i(0,0,2d)] 


(4-15) 


where  £  is  defined  in  (3-28) .  By  (2-42),  (4-3),  (4-11),  (4-12), 

(C-3),  and  (C-4),  we  have 


Z  -  Z 
NB,NB  1 1 


(4-16) 


In 


By  (2-42),  (4-3),  (4-10),  (4-t2),  (C-5),  and  (C-6) ,  we  obtain 


-jkL 

SNB,i  '  ~^-[jKL.0,rB)  -  4<I.,0,2d)  ] 


2  me 


{-Ci(ku  )+Ci(ku*  )+j[Si(ku,)-Si(ku’  )}) 


(4-17) 
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where  Ci  and  Si  are  defined  in  (C-7)  and  (C-9) ,  and 


u  -  L+  Jl2+i2 

2  *  B 


(4-13) 


u'2  -  L*  jL2+(2d)2 


(4-19) 


By  (2-42),  (4-1),  (4-11),  (4-12),  (C-5),  and  (C-6),  we  obtain 


Z  -  Z 
1  ,NB  MB, 1 


(4-20) 


Substituting  (4-2),  (4-10),  and  (4-12)  into  (2-42),  we  have 


{jkAl  [ij/(z  ♦•Lo.r  )  -  tfs  +-5,0,24)] 


ml  jue  m  —  m  2  B 


♦iKz  ♦|,0,r)-£(z  *|,0,2d) 

—  m+  2  B  —  m+  2 


-  iHz  +  ko,rn)  +  i|<(z  +  -5,0,24)}  m-2,3,...,SB-1  (4-21) 

—  n-  ^  B  m-  <- 


Substitution  of  (4-2),  (4-11),  and  (4-12)  into  (2-42)  gives 


.-4 


Z  MR  -  •2-rdjkA1  \Mz  -  5*0, r)  -if*  -  3*0,24) ] 
m,NB  jwe  m  —  m2  B  —  m2 


-4<V-2-0-‘b)  * 


*  ^Zn-'T’°'rB'  '  ^zn-”?,0’2d^  “"2.3 . n~'  <4-22> 


«•  V  -*  V  %■  V  V  V  '■*  */  *  •  */  ’  •  .  ’ ,  •  ,*•' ’/ 
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By  (2-37),  we  have 


Z  -  Z 

1m  ml 


m*2 ,3,«*» ,NB-1 


(4-23) 


Z  *  Z  m*2 ,3 , • • • ,  NB— 1 

NB,m  m,NB 


(4-24) 


The  evaluation  of  the  matrix  [z]  is  now  complete. 


A  A 

Finally,  we  evaluate  T  ,  T  ,  T  ,  and  T  for  m  ■  1,  2, 

ml  m,NB  1m  N3,m 

. . • ,  NA  as  follows. 


Substitution  of  (B-14)  with  z  "L/2  and  J“*J  into  (3-30)  gives 

2  ~  ~ 1 


^  2  /»c+  c+  L 

dl  e  x  z  ♦— 

T  ~  m _ J  m _  f  |  _  —  m  i  i  i  i  . 

ml  4tt  /  c+x2  2L  I7^I7°VV 

(xm  }  +d  J(xm  }  *d  +2> 


-jk  V  *c' 

_  m  m2  ra 

(x°-)2*d2 

m 


zc-*i 


/  C-x2  2  .  C-  lix< 

(xm  >  *d  *<z»  * 


,  -Jkj(x=-)^d'.(e=%i)2 
.  1  -  m  m2 

U0 


m  *  1j  2,  •  •  •  f  NA 


(4-25) 


where  $**,  x®*,  and  z®*  are  defined  in  (3-30)  and  (3-32). 

m  m  m 


Similarly,  substituting  (B-4)  with  z^-L/2  and  J^-J^  into  (3-30), 


we  obtain 

~jk— 

T  2  dlSL - I  — !!e _  f«. _ y~2 _ 


s* 


..  1/  c*\2  2  .  c*  Lv 2  -c- 

-jk4(x  )  +d  +(zffl  -  2)  xm 

e  +  — - - r - 5- 

<x«")  *d 


.C-  _  L 


/  c-»*  £.  ,  c-  t 

Pm  )  +d  +^zm  -  2 


-3k 


in*  If  2 )  •  •  •  f  NA 


(4-26) 


$  and  T  ,  for  m  *  1,2,...,NA,  can  be  obtained  by  (2-38). 

1m  NB,m 

By  solving  (2-14)  and  (2-15)  with  matrices  evaluated  above,  we 
can  obtain  the  current  distributions  in  the  aperture  and  on  the  wire. 


4.4.  Numerical  Results 


In  this  section,  numerical  results  for  the  magnetic  current  in 
the  aperture  and  the  electric  current  on  the  wire  are  presented.  To 
ensure  the  validity  of  our  formulation,  we  compare  our  numerical 
results  with  those  available  in  the  literature  for  several  examples. 
Very  good  agreement  is  obtained. 


The  aperture  of  the  first  example  is  an  electrically  small 
ellipse  with  major  axis  of  length  L  »0.02X  and  minor  axis  of  length 

BL 

W  *0.0021.  An  infinitely  long  wire  of  radius  r  -0.001X  is  located  at 

a  B 

y*d,  where  d*0.15X,  0.10X,  or  0.05X.  A  plane  wave  with  H*oa*  “e  ^Z/n 
ampere/meter  is  incident  from  region  a  of  Fig.  4.1.  Here,  n* /y/e.  We 
use  two  kinds  of  triangulations  to  model  the  aperture.  Figure  4.2(a) 
shows  a  triangulation  with  greater  patch  densities  at  the  extremities 
of  both  major  and  minor  axes  than  in  the  center  of  the  aperture. 


Figure  4.2(b)  shows  the  other  triangulation  which  has  greater  patch 

density  only  at  the  extremities  of  the  major  axis.  In  both  figures, 

NA,  NNODE,  NEDGE,  and  fJFACE  are  the  numbers  of  unkowns,  nodes,  edges, 

and  faces  of  the  aperture,  respectively.  About  the  triangulations, 

there  are  two  things  to  be  noticed.  First,  to  compensate  for  the  loss 

in  the  total  area,  we  should  put  the  boundary  nodes  outside  the 

aperture  so  that  we  get  the  correct  total  aperture  area.  Second,  to 

take  care  of  the  edge  effect,  we  need  a  higher  patch  density  around 

the  boundary  than  in  the  center.  In  Fig.  4.2,  the  boundary  nodes  are 

placed  on  the  boundary  of  an  ellipse  that  is  1 .01044  times  larger  than 

the  actual  one.  That  is,  21=1  .01044  T,  .  We  compare  our  results  from 

a 

the  two  triangulations  in  Fig.  4.2  to  illustrate  the  edge  effect. 
Table  4.1  shows  that  the  amplitudes  of  outward  traveling  TEF  currents 
on  the  wire  agree  very  well  with  those  obtained  by  using  the  formula 
in  f 1 ] .  It  also  shows  that  the  results  (NA*41)  using  the 
triangulation  of  Fig.  4.2(a)  have  better  agreement  than  those  (NA*17) 
using  that  of  Fig.  4.2(b). 

The  second  example  is  for  an  electrically  small  circular  aperture 

of  radius  r  *0.01A,  backed  by  an  infinitely  long  wire  of  radius 
A 

r  *0.001  A.  The  wire  is  at  y=d=0.03A,  0.04A,  0.05A,  or  0.06A.  The 

B 

excitation  is  a  plane  wave  with  H*oa«  -e”^kz/T1  ampere/meter.  Figure 
4.3  shows  the  triangulation  used  for  the  circular  aperture.  There  is 
a  higher  patch  density  around  the  boundary  than  in  the  center.  In 
order  to  get  the  correct  total  aperture  area,  we  use  i*i .0539rft  as  the 
radius  for  the  triangulation,  '"able  4.2  shows  that  the  amplitudes  of 
the  outward  traveling  TF.K  currents  agree  very  well  with  those  obtained 
by  using  the  formula  in  [ll.  For  examples  1  and  2,  our  results  also 


agree  with  those  of  the  formulas  in  l2J  and  |3'. 


The  third  example  is  for  a  narrow  slot  which  has  length  L  -0.5X 

a 

and  width  V  O.05\  and  is  centered  at  (0,0,0).  An  infinitely  long 
a 

wire  of  radius  r  =0.001A  is  located  at  y=d=0.1A.  A  plane  wave  is 
B 

i  08 

normally  incident  with  H  =1  ampere/ meter.  Figure  3  .7  shows  the 
triangulation  used  for  the  slot.  Fibres  4.  4,  4.5,  and  4.5, 

respectively,  show  that  the  total  x-directed  magnetic  current  M  in  the 
slot  and  the  real  and  imaginary  parts  of  the  electric  current  on  the 
wire  agree  very  well  with  those  in  f 4 .  figures  A. 5  and  4.5  show  that 
the  evanescent  current  on  the  wire  is  concentrated  in  the  region  for 
which  | z  j <0 . 25  X-  Beyond  this  region,  the  currents  are  two  outward 
traveling  TRW  waves.  Therefore,  the  choice  of  L=2.0A  used  in  our 
calculation  should  he  adeouate. 


Finally,  we  consider  a  circular  aperture  of  radius  r  “O.iA  bached 

A 

by  an  infinitely  long  wire.  The  wire  has  radius  r  ^O.OOU  and  is 

B 

located  at  y=»d=0.05A.  The  excitation  is  a  plane  wave  with 
H  *  -e  /ti  ampere/ meter.  V/e  use  the  triangulation  shown  in 
Fig.  4.3  for  the  aperture,  where  r=i  .Q53^r  is  the  radius  in  the 
patching.  Figure  4.7  shows  that  the  evanescent  current  on  the  wire 
drops  to  zero  rapidly,  and  can  be  neglected  beyond  a  small  region  of 
z*-0.25A  to  0.25A*  Note  that  both  the  evanescent  current  and  the 
outward  traveling  TRM  currents  are  discontinuous  at  z*0.  The  total 
current  is  continuous  everywhere  along  the  wire,  l’nis  implies  that 
pulses  must  begin  and  end  at  zm0.  The  replacement  of  delta  functions 
by  pulse  functions  in  Section  4.3  has  assured  continuity  of  the  total 
current  on  the  wire.  Figure  4.B  shows  the  radial  component  of  the 


magnetic  current  per  unit  length  V  crossing  a  circle  of  radius  5/6  r 

A 

in  the  aperture,  ''/hen  V  is  positive,  the  current  flows  away  from  the 
center  of  the  aperture.  When  V  is  negative,  the  current  flows  toward 
the  center.  It  is  seen  that  V  is  antisymmetric  about  the  z-axis 
(0*0.).  There  is  no  available  data  to  compare  with  for  this  example. 
In  Figs.  4.7  and  4.8,  H1Oa(0)  denotes  tfloa  at  z*0. 

X  x 
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Table  4.1. 


Amplitudes  I  and  I  of  outward  traveling  TEN  currents  on 

1  NB 

an  infinitely  long  wire  passing  by  a  small  elliptical 


aperture.  L  /X-0.02,  W  /X-0.002,  r/X-0.001,  and 

a  a  B 

ios  » 

”  “e  /n  ampere/metei.  (X-t  meter) 


Our  Solution 

Solution  of  [l]* 

NA  NB  L/X 

yl0“7 

w10'7 

ir/io"7 

W’0'7 

41  62  1.5 

-o.oogi-jo.1457 

-0. 0066- jO.  1303 

-j0.1554 

-jO.1474 

17  62  1.5 

-0.0066-j0.1391 

-0.0087-j0. 1  391 

41  62  1 .5 

-0.0077-j0.235B 

-0.0072-j0.2254 

-j0.2509 

-j0.2380 

17  62  1.5 

-0.0072-j0.2251 

-0. 0072- jO. 2251 

41  62  1.5 

-0.0078-j0.5401 

-0.0069-j0.5l62 

-jO.5763 

-j0.5476 

17  62  1.5 

-0.0072-j0.5153 

-0.0072- jO. 51 53 

1 

Nearly  the  same  as  the  solutions  by  the  formulas  in  [2]  and  [3] 
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f 

Fig.  4. 
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Fig.  4.6.  The  imaginary  part  of  the  current  on  an  infinitely  long 

wire  passing  by  a  slot  for  L  /x-0.5,  W  /X-0.05,  x  *0., 

a  a  c 

t  /X*0.001 ,  d/XO.1  ,  and  normal  plane  wave  incidence. 
B 
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Pig.  4.8.  The  radial  component  of  the  magnetic  current  per  unit 

length  crossing  a  circle  of  radius  5/6  r  in  a  circular 

A 

aperture  backed  by  an  infinitely  long  wire  with  plane  wave 


incidence.  r  /X-0.1,  d/X*0.05,  r  /X-0.001 ,  and 

A  B 


-9  *^Z/p  ampere/meter.  (M*56,  NB»62,  L/X*2.0f  and 


X  *1  meter) 
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Chapter  5 


APPLICATION  TO  THE  PROBLEM  OF  AN  ARBITRARILY  LOADED 
WIRE  BEHIND  AN  APERTURE  OF  ARBITRARY  SIZE  AND  SHAPE 


5.1 •  Problem  Specification 

In  this  chapter,  the  problem  is  generalized  to  that  of  an 

arbitrarily  loaded  wire  passing  by  an  aperture  with  excitations  coning 

from  either  one  or  both  sides  of  the  conducting  plane.  Figure  5.1 

shows  the  geometric  configuration  of  the  problem  to  be  considered.  An 

aperture-perforated  infinite  conducting  plane  of  zero  thickness  covers 

the  entire  x-z  plane  and  separates  regions  a  and  b.  In  region  a 

(y  <  0) ,  the  plane  wave  defined  in  (3-1)  is  incident  at  an  angle  (0  , 

o 

*  ).  In  region  b  (y  >  0) ,  a  z-directed  thin  wire  of  radius  r  is 
o  B 

terminated  by  loads  Z  and  Z  at  z  ■  -L  /2  and  L  /2,  respectively, 

LI  L2  w  w 

and  is  parallel  to  the  conducting  plane  at  a  distance  of  d.  There  are 
TEM  voltage  sources  v|  and  V*  applied  across  both  ends  of  the  wire. 
The  space  is  filled  with  a  loss-free  homogeneous  medium  of 
permeability  y  and  permittivity  e. 


The  generalization  is  to  include  the  reflection  effect  caused  by 
the  terminations  of  the  wire  into  the  moment  method  solution  developed 


in  Chapter  4.  The  objective  is  to  find  the  current  distributions  in 
the  aperture  and  on  the  wire.  We  then  evaluate  the  power  transmitted 
through  the  aperture.  In  addition,  an  equivalent  circuit  of  the 
aperture  for  the  transmission  line  mode  on  the  wire  is  obtained.  The 
derivation  of  such  an  equivalent  circuit  is  detailed  in  Appendix  D. 
Utilizing  the  circuit  and  transmission  line  equations,  currents  and 
voltages  of  the  TEM  mode  along  an  arbitrarily  loaded  wire  can  be 
calculated. 


5 .2.  Solution  Development 


As  in  Chapter  4,  the  presence  of  the  aperture  excites  two  outward 
traveling  TEM  currents  plus  evanescent  current  on  the  wire.  In 
addition,  there  are  reflections  of  these  currents  due  to  the 
terminations  of  the  wire.  The  wire  is  assumed  long  enough  that  the 


evanescent  currents  do  not  reach  the  terminations.  Thus,  the 
reflections  can  be  represented  by  two  pure  TEM  current  waves.  The 
situations  at  the  terminations  are  shown  in  Pig.  5.2,  where 
(I1®Jk\,INBe“JkZuz)  are  the  outward 

traveling  TEM  currents  and  the  pure  TEM  currents,  respectively,  and 
the  arrow  denotes  the  direction  of  the  propagation.  The  original 
problem  in  Pig.  5.t  is  then  equivalent  to  that  in  Fig.  5*3,  where  the 

loads  and  the  voltage  sources  are  replaced  by  the  pure  TEM  currents 
♦  “jfcz  .  “  jkz 

_  ®  ^2  an(1  _  *1  e  Uj5  traveling  on  an  infinitely  long  wire. 
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ft 

ft 


Here,  I  and  I  are  related  to  the  voltage  sources,  the  loads,  and  the 
outward  traveling  TEM  currents.  The  relationships  are  derived  as 
follows. 

Referring  to  Fig.  5.2,  the  voltages  at  the  terminations  are 


-  vj  -  (5-° 


♦  T  J*z\ 


7  (t*  t"  jkz  r  “jkz\  „i  a.  „  fr*  -jk®  ~  jkz  -jkzx  . 

V1  e  -I  e  *lHBe  )  -  V2  ♦  yi  e  +1  eJ  *Imm  )  (5- 


LNB" 


-2) 


where  Z^  is  the  characteristic  impedance  of  the  transmission  line 
formed  by  the  wire  and  the  ground  plane. 


At  z-0,  (5-1)  and  (5-2)  become 


VI+  +  r  +  ii> 


VI+  -  I_  +  W  ’  V2  *  Z2(I+  +  I_  +  W 


(5-3) 


(5-4) 


Here,  and  are  V*  and  V*  referred  to  z"0,  respectively.  Z^  and 

Z  are  the  impedances  Z  and  Z  referred  to  z*0,  respectively. 
c.  LT  L2 


By  solving  (5-3)  and  (5-4)  for  I*  and  l“,  we  obtain 


I+  *  To  *  Vi  *  VHB 


i"  -  iZ  *  c0i.  ♦  c  i 

0  21  3  NB 


(5-5) 


(5-6) 
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where 


(5-7) 

(5-8) 

(5-9) 

(5 -VO) 

(5-11) 

(5-12) 

(5-13) 


and  r2  ere  the  reflection  coefficients  at  z*0”  an(j  o+, 

respectively.  Note,  I*  and  i"  are  infinite  only  when  Z^m 00  or 

Z^-Z^.  The  first  ease  (Z^Z^*00)  does  not  exist  for  the  loaded  wire. 

The  second  case  (Z  »-Z  )  occurs  when  the  resistance  of  one  of  the 

1  2 

loads  is  negative  or  the  resistances  of  both  loads  are  zero.  However, 
usually  the  loads  are  passive  and  the  sources  are  lossy.  Thus,  this 
caae  does  not  arise.  Therefore,  I+  and  I~  are  finite  . 
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By  (5-5)  and  (5-6),  we  have 


a 


-  4}  +  X1-1  +  Imr*B 

I  z  |  <  00 

(5-14) 

-  "  +  xi-i  +  tnAb 

1  Z  {  <  oo 

(5-15) 

where 

4  ■  v'jk4 

|  Z  !  <  oo 

(5-16) 

4  ■  VJk\ 

|  Z  |  K  oo 

(5-17) 

4  •  VJk\ 

J  Z  )  <  0° 

(5-18) 

5  ■  VJk  % 

|  z|<  00 

(5-19) 

4  •  V3k4 

|z|<  00 

(5-20) 

4b  -  v*’ % 

|  Z  1  <  oo 

(5-21) 

We  now  consider  the  problem  in 

Pig.  5.5. 

There  are  an  equivalent 

magnetic  current  sheet  M  over  the  aperture,  the  pure  TEM  currents,  and 
an  induced  electric  current  £  on  the  wire  due  to  M.  We  can  use  the 
same  expansion  functions  for  M  and  £  as  those  in  Chapter  4.  The 
formulations  and  matrices  developed  in  Chapter  4  can  also  be  utilised 
with  some  modifications. 
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tv 


s, 
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We  now  apply  the  boundary  conditions  stated  in  Section  2.2  to  the 
problem.  The  first  condition  that  the  tangential  _E  is  continuous  at 
the  aperture  is  assured  by  the  magnetic  current  sheets  jl  and  -JI  over 
the  two  sides  of  the  aperture  region.  Due  to  the  additional  pure  TEM 
currents,  the  second  condition,  (2-9),  becomes 


+1  J+  +1  J"  ) 

1-»  NB-NB  NB-NB' 


on  A 


(5-22) 


Here,  J  consists  of  two  outward  traveling  TEM  currents  plus  evanescent 
current.  The  condition  (2-10)  remains  the  same  because  the  tangential 
J2  from  a  pure  TEM  current  vanishes  on  the  wire. 


Therefore,  the  matrices  [ya+Yb],  [Z],  fT],  tia,  and  Vib  remain 

the  same  as  those  in  Chapter  4  while  [t]  and  tib  muat  be  modified. 
Note  that  here  [t]  is  defined  in  (2,-*8)  or  (2-40)  with  replaced  by 

and  4b  by  4b^b+5b  ’  and  ?ib  i3  defined  aa 


(5-23) 


We  now  modify  the  [t]  matrix.  By  (2-40)  and  (5-22),  the  modified 
matrix  is  the  sum  of  the  [t]  due  to  J  (same  [t]  developed  in  Chapter 
4)  plus  additional  terms  T®,  and  T®>NB  due  to  J*.  J_~, 
for  m*1 ,2, • . . ,NA.  The  additional  terms  are 


4b-  and 


’St  • 


(5-24) 


(5-25) 


where  J_*  *  and  are  the  images  of  and  J ,  respectively. 


By  [4,  Eqs.  (80),  (83)]  and  (3-29),  the  tangential  magnetic 

♦  lies 

field  in  the  aperture  due  to  a  pure  TEM  current  e  **  u  on  the  wire  is 

— z 


Hr  -  — %  Tu 

"°  2>(x2*42f* 


(5-26) 


By  (5-t7),  (5-18),  (5-20),  (5-21),  and  (5-26),  we  have 


V4'4-i'U'7>  •  2<°,^  *  c£> 


(5-27) 


n-i  m>  *  2(cb4  * 


(5-2B) 


Substitution  of  (3-4),  (5-27),  and  (5-28)  into  (5-24)  and  (5-25)  gives 


Ta.  «  C  T’  ♦  CT+  m  -  1,  2 . HA 

mi  1  m  2  m 


(5-29) 


t:,wb  -  v:  ♦  c3t;  . na 


(5-50) 


where 


tjkzj'  tjkz0- 

di  xc+e  m  ;c%  m 


.  c+»2  2  ,  c-»2  2 

<xm  )  +d  (xnr  >  +d 


are  defined  in  (3-30)  and  (3-32) 
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By  (3-4),  (5-16),  (5-19),  and  (5-23),  the  mth  element  of  the 

V  +  . 

source  vector  I1D  due  to  io  and  Jq  is 

llb  "  IITI  +  Ol  m-1,2 . NA  (5-32) 

m  0  m  urn 


The  evaluation  of  the  matrices  is  now  complete*  By  solving 
(2-14)  and  (2-15)  with  the  matrices  evaluated  above,  the  currents  in 
the  aperture  and  on  the  loaded  wire  can  he  obtained. 


The  total  current  on  the  wire  is 


J+  J  ♦  J 


Thus,  the  total  TEH  currents  at  z*0  and  0  are 


ITEH 

ITEM 


(l-vyi,  +  (VC3)rNB  *  T0  +  T0 
+  (lfVC3)XNB  +  V  ^ 


(5-33) 


(5-34) 

(5-35) 


It  is  noteworthy  to  verify  that  the  problem  considered  in  Chapter 
4  is  a  special  case  here.  The  infinitely  long  wire  considered  in 
Chapter  4  is  equivalent  to  a  loaded  wire  with  Zi*z2*Zo  and  •(). 

Thus,  (5-5)  -  (5-13)  give  I  -1  «o.  The  current  on  the  wire  is  then 

reduced  to  that  of  Chapter  4. 
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Finally,  for  a  loaded  wire  whose  length  is  very  short,  the 
solution  developed  here  may  not  apply •  For  this  case,  the  method  of 
Chapter  3  can  be  utilized.  However,  the  pulse  functions  on  the  wife 
must  be  extended  to  the  terminations.  In  addition,  the  impedance 
matrix  [z]  must  be  modified  by  adding  a  diagonal  matrix  with  loads  as 
its  diagonal  elements  corresponding  to  the  terminations. 


3«3«  Power  Transmitted  through  the  Aperture 

In  this  section,  we  derive  the  power  transmitted  through  the 
aperture  into  region  a  of  Fig.  5.1  when  the  excitation  is  incident 
from  region  b. 

By  [7,  Eq.(l-57)],  the  complex  power  transmitted  through  the 
aperture  into  region  a  is 
P®  -  -  || (EaxHa*)-UydS 

-  -  |J(uyxEa).Ha*dS 


y 


NA 

NA 

_  *  f  f 

l 

l  V  v  1 

L  m  n  )  j  • 

m*1 

n*1  A 

NA 

NA 

I 

l  v  vV* 

ro  n  inn 

m*1  n*1 


70 


1^1 


*[Ya]V 


(5-36) 


where  "  *  M  denotes  complex  conjugate.  In  (5-36) ,  the  vector  identity 

(AxB)  •  £  *  (£xA)»3,  (2-1)  with  n-u_^,  (2-11),  and  (2-t6)  are  used.  The 

time-average  power  P  transmitted  through  the  aperture  is  given  by  the 

t 

real  part  of  Pa. 

f 


3.4.  Evaluation  of  an  Equivalent  Circuit 


In  this  section,  an  equivalent  circuit  at  2*0  of  the  aperture  for 
the  transmission  line  mode  on  the  wire  is  presented.  For  simplicity, 
a  plane  wave  incident  from  region  a  of  Fig.  5.1  is  the  only 
excitation.  Using  this  circuit,  we  can  obtain  the  TEM  currents  and 
voltages  on  the  wire. 


i 


m 


The  derivation  of  the  equivalent  circuit  is  detailed  in  Appendix 
D.  In  Appendix  D,  the  solution  for  the  current  on  a  wire  with  matched 
loads  (or  infinitely  long  wire)  is  utilised  to  evaluate  the  impedances 
and  sources  of  the  circuit.  Ve  therefore  consider  the  wire  as 
infinitely  long. 


The  general  T  network  is  shown  in  Fig.  D.5.  The  impedances  Z^ , 

Z®,  and  Z®  defined  in  (D-21)  -  (D-23)  depend  on  Z  ,  I  ,  I  ,  I  ,  and 
c  v  0  1  NB  1 

*  ^  ^ 

I  .  Here,  ( f  ,1  )  and  (I  fl„n)  are  the  amplitudes  of  outward 
traveling  TEH  currents  in  the  -z  and  the  *z  directions  on  the  wire 
when  the  original  incidence  is  replaced  by  TEM  current  excitations 
and  e*^®^  applied  to  the  wire,  respectively.  The  sources  V® 
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and  Vj  in  (D-24)  and  (D-25)  depend  on  ZQ,  Z®,  Z®,  Z®,  I, ,  and  1^.  Iy 

and  I  are  the  amplitudes  of  outward  traveling  TEM  currents  in  the  -z 
NB 

and  the  +z  directions  on  the  wire,  respectively,  when  the  excitation 
ia  the  original  incident  field. 


If  the  aperture  is  symmetric  about  z*0  or  is  small,  the  network 
in  Pig.  D.7  is  used.  The  impedances  Z®  and  Z®  and  sources  V®  and  V® 

A 

defined  in  (D-28)  -  (D-31)  depend  on  Z  ,  I  ,  I  ,  I  ,  and  I„_. 

0  t  NB  1  NB 

If  the  aperture  induces  symmetric  outward  traveling  TEM  currents 

on  the  wire,  the  network  in  Fig.  D.8  with  impedance  Ze  in  (D-32)  and 

6  A 
source  V  in  (D-33)  is  utilized.  The  network  depends  on  ZQ,  1^,  and 

V 

If  the  aperture  induces  antisymmetric  outward  traveling  TEM 
currents  on  the  wire,  the  network  in  Fig.  D.9  with  impedance  Za  in 
(D-24)  and  source  V®  in  (D-35)  is  used.  The  network  depends  on  Zq» 

A 

It,  and  1^ . 

To  obtain  the  network,  we  first  calculate  the  amplitudes  of  the 
outward  traveling  TEM  currents  as  follows. 

/V  A 

To  calculate  (1^,1^),  (^,1^),  and  (I1*INB)t  solve  (2-14) 

and  (2-t5)  with  the  matrices  [Ya+yb],  [z],  [t],  and  [f]  evaluated  in 
Chapter  4  and  with  for  the  current  on  the  wire.  For  (l>,t^g), 

the  source  vectors  Ii8*0  and  Iib  defined  in  (5-32)  with  1^*1  and  1^-0 
are  used.  For  ( 1^ ^ jjg) *  tia«0  and  I defined  in  (5-32)  with  !*■()  and 
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I““1  are  used.  For  (I  ,1  ),  I*a  defined  in  (3-35)  and  1*^*0  are 

U  *  NB 

used. 

Next,  we  substitute  these  amplitudes  and  into  the  expressions 
for  the  sources  and  impedances  of  the  network. 


5*5.  Numerical  Results 

Numerical  results  for  the  current  distributions  in  the  aperture 
and  on  the  wire  are  presented  in  this  section.  Numerical  results  for 
the  power  transmitted  through  the  aperture  are  also  given  here.  In 
addition,  we  calculate  the  impedances  and  sources  of  the  equivalent 
circuit.  The  TEM  current  on  the  wire  is  then  computed  from  this 
circuit  and  compared  with  that  obtained  from  the  solution  of  Section 
5.2. 
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Next,  we  consider  a  slot  of  length  L  «0.5X  and  width  W  "0.05X 

a  a 

backed  by  a  wire  of  radius  r  «0.001X.  The  wire  is  located  at  various 

6 

positions  from  y“d*0.1\  to  0.4X.  A  plane  wave  is  normally  incident 

from  region  a  of  Pig.  5.1.  Pigure  3.7  with  x  *0.  is  used  for  the 

c 

triangulation  of  the  slot.  Figures  5.5  and  5.6  show  the  source  and 

impedance  of  the  equivalent  circuit  (in  Pig.  D.ft)  of  the  transmission 

line  mode  on  the  wire.  Our  solution  agrees  very  well  with  that  in 

[4].  From  this  circuit,  we  can  also  calculate  the  TEM  current  at  z*0 

on  the  wire,  by  terminating  the  two  ports  of  the  circuit  with  Z^  and 

Z^.  This  is  shown  in  Fig.  5.7.  The  TEN  current  at  z“0  is  equal  to 

the  loop  current  I  in  Fig.  5*7.  Table  5.1  gives  the  results  for 

d*0.1X  and  various  loads.  It  shows  that  the  results  are  (almost) 

identical  to  those  obtained  by  the  solution  in  Section  5.2.  They  also 

agree  well  with  those  in  [6].  In  addition,  the  results  of  the  last 

three  cases  show  that  the  total  TEM  current  at  zm0  is  unchanged  if  the 

sum  of  loads  Z  and  f»  remains  unchanged. 

1  2 


.*• 


& 


200 


XA  OUR  SOLUTION  (NA*I9) 
-  SOLUTION  IN  14] 


Fig.  5*5*  Equivalent  source  for  the  TEM  mode  at  2*0  on  an  arbitrarily 
loaded  wire  passing  by  a  slot  with  normal  plane  wave 
incidence.  L  /X*0.5»  W  /X  *0.05,  r_M  *0.001,  and  Hioa-1 

a  a  B  * 

ampere/meter.  meter) 


Table  5«1.  Total  TEM  current  at  z-0  on  a  loaded  wire  passing  by  a 

slot  with  normal  plane  wave  incidence.  L  /X  - 0.5, 

a 

W  /X -0.05,  d/X-0.1,  x  /X -0.001,  and  Hioa«1 

&  D  X 

ampere/meter.  (  X -1  meter,  NA-19,  NB-62,  and  L^/X-1.5) 


Reflections 

Solution  of 
(5-34)  or 

Calculated 

Solution 

F1 

r2 

(5-35) 

from  Pig.  5.7 

in  [6] 

0.1 

I. 

0 

0.2191 -j0.0792 

0. 2192- jO. 0792 

0.2062- jO. 0861 

1.0 

0 

0 

0 

0 

B 

0 

0.02R7-j0.01 29 

0.0287- jO.OI 29 

0.0264-j0.0136 

0.8182 

0.8182 

0.0287-j0.0129 

0.0287-J0.0129 

0.0264-j0.0136 

0.8750 

0.6667 

0.0287-J0.0129 

0. 0287-jO. 01 29 

0.0264-j0.0136 
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Chapter  6 


CONCLUSION 

A  moment  method  solution  for  the  problem  of  electromagnetic 
coupling  between  an  aperture  and  a  conducting  body  is  developed.  Both 
the  aperture  and  the  conducting  body  are  arbitrarily  sized  and  shaped. 
This  method  is  then  implemented  in  specific  cases  whereby  the 
conducting  body  is  a  wire.  The  wire  is  finite-length  unloaded, 
infinitely  long,  or  arbitrarily  loaded.  The  aperture  is  modeled  by 
triangular  patches,  which  can  conform  closely  to  any  geometry  and 
permit  greater  patch  densities  on  those  portions  where  more  resolution 
is  desired.  Local  position  vectors  are  used  as  the  expansion 
functions  in  the  aperture.  There  are  two  outward  traveling  TEM 
currents  and  evanescent  current  on  the  infinitely  long  wire.  The 

evanescent  current  is  found  to  exist  in  a  small  finite  region  near  the 
aperture.  Numerical  results  show  very  good  agreement  with  available 
data  in  the  literature. 

Further  work  is  recommended  for  more  general  proolems  as  follows. 
First,  consider  the  problem  of  a  set  of  arbitrary  apertures  with  a  set 
of  parallel  wires.  Each  wire  can  be  finite-length  unloaded, 
infinitely  long,  or  arbitrarily  loaded  and  can  be  in  different  media. 
One  can  use  the  came  expansion  functions  as  those  developed  in  this 
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work  for  each  aperture  and  each  wire.  The  matrices  in  (2-14)  and 
(2-1 5)  must  be  modified  to  include  the  coupling  among  apertures  and 
wires.  In  addition,  different  permeabilities  and  permittivities  must 
be  included  in  the  matrices.  By  solving  (2-14)  and  (2-15),  the 
currents  in  the  apertures  and  on  the  wires  can  be  obtained. 

The  problem  of  a  set  of  arbitrary  apertures  with  a  set  of 
arbitrary  conducting  bodies  can  also  be  done.  To  solve  this  problem, 
one  can  use  the  formulation  developed  in  Chapter  2  and  the  matrix 
equations  (2-14)  and  (2-15).  Triangular  patching  is  recommended  to 
model  apertures  and  conducting  bodies.  Local  position  vectors  can  be 
used  as  expansion  functions  .  The  work  may  be  complicated,  but  it  can 
be  done  without  too  much  difficulty. 
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APPENDIX  A 


SURFACE  INTEGRATIONS  FOR  FIELD  PROBLEM 


In  this  appendix,  we  evaluate  the  following  surface  integrals: 


1  1-U 


■  I  ! 

n-o  5«o 


-jkn 

e—dUri 


i  i-n 


■t*  I  I 

n-o  £-o 


(A-2) 


i  i-n 


V  J  I  ’V** 

n*o  c*o 


(A-3) 
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Here,  R-|r^-r' |  is  the  distance  between  a  source  point  at  _r'  in 
triangle  T  and  an  observation  point  at  jr.  The  area  coordinates  of 
triangle  T  are  £  and  n  defined  in  (3-iO  and  (3-*2). 


To  avoid  the  singularity  of  these  integrals  at  _r-:r' * ,  we  rewrite 
(A-l)  -  (A-3)  as 


1  1-r 
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i>0  £-0 
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The  first  terms  of  (A-4)  -  (A-6)  can  be  approximated  by  using  a 


seven-point  Gaussian  quadrature  formula  developed  for  a  triangular 


region  [l-6],  The  formula  is 


t  irn 

j  I  e(l*5»n  )<*£<*) 

0-0  k-  0 


i  Vt*(i,V'i) 

i«i 


(a-7) 


Here,  g  is  the  integrand,  (5  »n^)  is  a  point  in  triangle  T,  is  the 


weight,  and  a^  is  the  coefficient.  They  are 


wi-0.5 


i-1  ,2 . 7 


a^-0.225 


a  -a  -a  -0.1323942 
2  3  4 


a  -a  -a  -0.1259392 
5  6  7 


5,-^-0.333333 
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£  .£  «n  -n  -0.470U206 
3  4  2  4 


£  -ri  -0.7974269 
5  6 


£  •£  — ri  -ri  -0.1012865 
^6  7  5  7 


To  evaluate  the  second  terms  of  (A-4)  -  (A-6),  we  first  introduce 


some  vectors  and  variables  as  shown  in  Fig.  A.1.  0  denotes  the  origin 


of  the  global  coordinates.  £  is  the  projection  of  _r' onto  the  plane 


of  the  triangle  T.  The  distance  between  the  observation  point  and  the 


plane  of  the  triangle  T  is  d.  The  contour  along  the  boundary  of  T 


is  defined  as  C,  which  is  the  sum  of  three  straight  lines  ,  C^,  and 


C  along  the  edges  of  the  triangle,  u  ,  u  ,  and  u  are  the  outward 
3  “2  “3 


unit  vectors  normal  to  each  edge.  Point  P  denotes  the  projection  of 


the  observation  point  onto  T.  For  edge  i  (i-1,2,3),  P  and  P  are 

~i  — oi 


the  vectors  from  point  P  to  a  point  on  and  to  the  projection  of  the 


observation  point  onto  line  C  ,  respectively.  The  unit  vectors  of  P 

i  ~i 


and  o  are  u  and  u  .  Point  9  denotes  the  projection  of  point  0 
•“oi  ~Pi  -Poi 


onto  line  C^.  Variable  1^  is  the  length  measured  from  point  Q  to  a 


point  on  and  is  in  the  direction  of  a  unit  vector  _u  .  This  unit 


vector  is  in  the  direction  of  the  boundary  contour  of  triangle  T  along 


edge  C^.  The  end  points  of  and  the  projection  of  the  observation 


point  onto  C  are  at  1  -l”,  1.,  and  1  . ,  respectively, 
i  i  i  i  oi 


Ve  now  evaluate  the  second  term  of  (A-4),  which  can  be  written  by 


using  (3-16)  as 


1  2A 


Hk¥t 


(A-9) 


.■  >  >  -ji  "  i 


V’  *f(p) 

S  — 


(A-12) 


where  u  is  the  unit  vector  in  the  p  direction,  and  V'  is  the 

M  Q 

surface  divergence  with  respect  to  the  coordinates  of  jr' . 

By  (A-ll)  and  (A-12),  we  can  evaluate  (A-9)  using  the  formula  in 
(A-IO).  However,  note  that  application  of  (A-10)  requires  the 
continuous  differentiability  of  the  integrand  on  the  domain  of 
integration,  S.  This  is  not  the  case  when  P  is  in  S  or  on  C. 
Therefore,  a  region  T£,  defined  as  the  intersection  with  T  of  an 
infinitesimal  sphere  with  radius  c  centered  at  P,  is  isolated  for 
separate  treatment.  The  boundary  of  T^  j.s  defined  as  C^,  anda(r.)  is 
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the  angular  extent  of  C£  about  P.  The  region  of  T  excluding  T£  is 

T-T  ,  and  the  corresponding  contour  is  defined  as  C-C~.  Thus,  (A-9) 

£  ^ 

is  rewritten  as 


The  second  integral  in  (A-*3)  is 


lim  cx(r)C>je2+d^-  |d  |)  -  0 

£-*0 


(A- *3) 


(A-14) 


By  (A-1‘0)  -  (A-14),  we  have 
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Here,  a(0  can  be  written  as 
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where 
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-  Jd2+(poi)2+uf)2 


(A-23) 


Substitution  of  (A-16)  and  (A-22)  into  (A-15)  gives 
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?A  A  — Oni  Pnl 


2A  1  ■=poi 
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p  ai: 
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A  +Poi+U|Ri 
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In  (A-24), 


tan(A-B) 


tanA-tanB 
1*tanA  tanB 


(A-25) 


is  used. 


There  are  four  things  to  be  noted  about  (A-24);  (1)  when  p  «0, 

01 

♦  + 

the  In  term  vanishes,  (2)  when  p  .^0,  the  arguments  R“+A1“  of  the  In 

oi  i  i 

term  are  positive,  (3)  when  d-0,  the  tan"1  terms  vanish,  and  (4)  when 

1  2  2  + 

dj«0,  the  arguments  of  tan"  are  finite  (i.e.,  d  +Poi+ldlRi  >  °)  •  In 
addition,  (A-24)  can  be  applied  to  any  surface  with  multiple  edges. 
If  there  are  N  edges,  the  summation  index  i  is  up  to  N. 


Finally,  we  consider  the  second  terms  of  (A-5)  and  ( A— 6) •  They 
have  been  evaluated  in  [l2,  Appendix  c]  in  terms  of  the  value  of  1^ 
defined  in  (A-9)  by  using  [l9,  Eqs.  (380.011),  (380.201)].  The 

results  are 
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Here,  r  ,  r  ,  and  r  are  defined  in  Fig.  3.6. 

— 1  ~2  ~3 

The  evaluations  of  the  integrals  in  (A-l)  -  (A-3)  are  now 

complete. 


OBSERVATION 

POINT 


Geometrical  quantities  associated  with  an  observation 
point,  a  source  point,  and  edges  of  source  triangle  T. 
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APPENDIX  B 


ELECTRIC  AND  MAGNETIC  PI ELDS  DUE  TO  OUTWARD  TRAVELING  TEM  CURRENTS 


Here,  we  derive  fields  due  to  two  outward  traveling  TEM  currents 
starting  from  z*z  and  -z  and  traveling  to  infinity  in  the  z  and  the 

r  2 

-z  directions,  respectively.  The  currents  are  defined  as 
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We  first  evaluate  the  ♦-component  of  the  magnetic  field  due  to  J* 
at  a  point  (p,  <j>,z)  as  follows. 
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Here, 
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(B-5)  is  used  on  the  third  step,  and  the  relationship 
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are  used  on  the  last  step. 


Similarly,  the  ^-component  of  the  magnetic  field  due  to  J“  can 
be  obtained  as  follows. 
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Comparing  (B-11)  with  (B-3)  and  (B-13)  with  (B-5),  we  can  obtain 
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The  z-conponents  of  electric  fields  due  to  J+  and  J*  are  then 
obtained  by  the  equation 
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(B-16) 


Substituting  (B-4)  and  (B-14),  respectively,  into  (B-16),  we  obtain 
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APPENDIX  C 


DERIVATIONS  OP  INTEGRAL  FORMULATIONS  USED  IN  CHAPTER  4 

The  integral  formulations  used  in  Chapter  4  are  listed  as  follows 
and  proved  later  in  this  appendix* 
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P  is  a  pulse  function  over  a  very  small  region  A l  and  P>  0. 


We  now  start  to  prove  (C-t).  Letting 
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,  we  have 


du  •  k(  )ds 

4z  *p 


(C-1?) 


Substitution  of  (C-f2)  and  ( C — 1 3 )  into  the  left-hand  side  of  (C-1) 
gives 


00  00  00 

1 4“‘“-  \  J  **?*■-  1 2j^d“  i 


-Ci(kp)  -  j[L-  Si(kp)] 


Thus,  (C-1)  is  proved. 


vV.-VvVvvvVvv'v 
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The  left-hand  side  of  (C-2)  can  be  rewritten  as 


By  (C-12)  and  (C-13),  we  have 


-  -  Ci(kUj )  -  j[J--  Si(kUi ) ]  (C-15) 

Using  (C-l),  (C-14),  and  (0-13),  we  obtain  the  right-hand  side  of 
(C-2). 

Replacing  z  by  z'+a  in  the  left-hand  side  of  (C-3)  and  using 
(C-t),  we  obtain  the  right-hand  side  of  (C-3). 


The  left-hand  side  of  (C— 4)  is  equal  to 


The  first  term  equals 
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-4irAl  e"^ajj>(0#0#P) 

Approximating  the  integral  with  respect  to  z*  in  the  second  term  by 
the  product  of  Al^  and  the  integrand  at  z'"0,  and  then  using  (C-3),  we 
have  the  second  term  equal  to 

jkAloe-jka|-Ci(kP)-j[y-Si(kP)]{ 

We  thus  complete  the  derivation  of  (C-4). 

Replacing  z  by  z'-a  in  the  left-hand  side  of  (C-5)  and  using 
(C-l)  and  (C-2)  with  a  replaced  by  2a,  we  obtain  the  right-hand  side 
of  (C-5). 

Finally,  using  the  same  procedures  as  those  in  deriving  (C-4)  and 
(C-5),  we  obtain  (C-6). 


APPENDIX  D 


EQUIVALENT  CIRCUIT  OF  THE  APERTURE 
FOR  THE  TRANSMISSION  LINE  MODE  ON  A  LOADED  WIRE 

In  this  appendix,  we  derive  formulations  for  an  equivalent 
circuit  of  the  aperture  for  the  transmission  line  mode  at  z*0  on  a 
wire.  The  wire  is  in  the  z-direction,  terminated  with  arbitrary 
loads,  and  parallel  to  an  aperture-perforated  infinite  conducting 
plane.  The  aperture  is  arbitrarily  sized  and  shaped  and  is  centered 
at  z*0.  Incident  fields  come  from  either  one  or  both  sides  of  the 
conducting  plane.  This  circuit  can  be  utilized  to  calculate  the  TEM 
currents  and  voltages  on  the  wire. 

It  is  known  that  an  equivalent  circuit  of  a  wire  passing  by  an 
aperture  is  a  two-port  network  as  shown  in  Fig.  D. i.  The  parameters 

[z°]  depend  only  on  the  geometry  of  the  problem.  They  can  be  obtained 
by  removing  the  incident  fields  from  the  system  and  applying 
mathematically  arbitrary  excitations  to  the  transmission  line  formed 
by  the  wire  and  the  conducting  plane.  The  sources  V®  and  V®  can  be 
obtained  from  the  incident  fields  and  [z®].  in  order  to  find  the 
equivalent  circuit,  the  solution  for  the  current  on  an  infinitely  long 
wire  (or  a  matched  wire  with  the  characteristic  impedance  Z^)  is 
needed.  Therefore,  we  now  consider  the  wire  as  infinitely  long. 


(A)  Network  Parameters  [Ze] 


I* 


l* 


To  obtain  the  elements  of  [ze],  we  use  the  definitions 


v 

1 


♦  Z 


12 


i 

2 


(&-1) 


v 

2 


21 


♦  Ze  i 
22  X2 


(D-2) 


where  i^ ,  i^,  ,  and  are  port  currents  and  voltages  at  port  1 

(s*0  )  and  port  2  (z*0  ),  respectively,  with  reference  directions 

shown  in  Fig.  D.2. 


To  represent  an  excitation,  we  mathematically  apply  a  TEM  current 

wave  I  e  *^Zu  to  the  transmission  line.  I  la  an  arbitrary  constant. 
— z 

This  excitation  will  induce  an  equivalent  magnetic  current  sheet  in 

the  aperture.  The  magnetic  current  sheet  then  excites  two  outward 

traveling  TEM  currents  I  eJ  u  and  I  e  J  u  traveling  on  both 

1  "T!  NB  ~Z 

semi-infinite  halves  of  the  wire  .  It  also  excites  evanescent  current 
on  the  wire  in  the  vicinity  of  the  aperture.  Figure  D.3  shows  the 
propagation  directions  of  the  excitation  and  the  outward  traveling  TEM 
currents. 

Combining  Figs.  D.2  and  T>.3,  we  obtain  port  currents  i^  and  i^ 
and  port  voltages  ^  and  v^.  They  are 


A 


♦ 

I 


A 


♦  I, 


(D-3) 
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A 


(D-4) 


V1  -v 


*2  *  V1*  *  *1B> 


(D-5) 


(D-6) 


Here,  Z  «601n(2d/r  )  is  the  characteristic  impedance  of  the 
0  B 

transmission  line,  d  is  the  distance  of  the  wire  to  the  conducting 
plane,  and  r^  is  the  radius  of  the  wire.  Since  [ze]  i3  independent  of 
the  excitation,  (D-1)  and  (D-2)  hold  for  the  port  currents  and 
voltages  defined  in  (D-3)  -  (D-6).  That  is, 


/s 


Zuh  *  Z12% 


(D-7) 


* 


Z®  t 
21  1 


(D-8) 


To  solve  for  [Ze],  two  more  equations  relating  the  port  currents 
and  voltages  to  [ze]  are  needed,  in  addition  to  (D-7)  and  (D-8).  One 
can  apply  a  TGM  current  excitation  to  the  wire.  Again,  the 
aperture  excites  two  outward  traveling  currents  1^ e°  ^  and 
I^B*  J  Jig  traveling  on  both  semi -infinite  halves  of  the  wire,  as 
shown  in  Fig.  D.4.  The  port  currents  and  voltages  for  this  case  can 
be  obtained  by  combining  Figs.  D.2  and  D.4,  giving 


I  ♦  I< 


(0-9) 


1JWLWII U  L  « 'J n  W  W- I?  l5J 


* I >;  ^Xv fTf>  r* rv»V>\^ > 


106 


A  ^ 

f2  •  -(1‘  *  *»> 

(D-10) 

•  -vr  ♦  I*,> 

(D-U) 

/\ 

^  A 

v  -  Z  (I  -  I  ) 

2  0V  NB  ; 

(D-12) 

Similar  to  (D-7)  and  (D-8),  we  have  a  pair  of  equations 

port  currents  and  voltages  to  [ze],  They  are 

relating  these 

v  «  7®  i  *  7®  i 

1  11*1  12*2 

(D-»3) 

1?  «  Z®  i  ♦  Z®  i 

2  21  1  22  2 

(D-14) 

For  simplicity,  I  and  I  are  chosen  to  be  the  constant  t.  Thus, 

(D-3)  -  (D-14)  give 

z?>  ■  [(,*tra)('*V‘(,-i,)('*v]p 

(D-15) 

Z*2  -  2(l*!t)F 

(D-16) 

z|,  ■  z('-iH)r 

(D-17) 

z22  ■ 

(D-18) 

where 
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F 


z 

_a 


1  )-(i+f  )(ui  ) 
NB  NB  1 


(D-19) 


By  the  reciprocity  theorem  [7,  Sec.  3-8],  we  have 


(D-20) 


The  network  is  then  equivalent  to  a  T  network  as  shown  in 
Fig.  D.5,  where  the  impedances  are,  as  obtained  by  using  ( D— 1 5)  - 
(D-20), 

z*  *  z?rz?2  “  f(l4i)(uV’(1‘!NB)(,^i)]F  (D“20 

Z2  *  Z*2  ’  Z21  "  2(l4NB)F  (D_22) 

Z3  *  Z22'Z2J."  t(l41)(l4NBMUiNB)(l"VlP  (D_2,) 


(B)  Sources  V*  and 

If  both  ports  of  Fig.  D.5  are  terminated  with  matched  loads,  as 
shown  in  Fig.  D.6,  the  network  corresponds  to  our  original  problem. 
The  original  incident  fields  induce  two  outward  traveling  TEM  currents 
I  e^Bu  and  I  e”^*  on  the  wire.  The  loop  currents  i  and  i  in 
Fig.  D.6  are  equal  to  It  and  respectively.  Thus,  Fig.  D.6  gives 


maun  w.iwim.t'.-  .*•  • 


7zl^v~, 


v» '  h*i’  •  : 
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the  sources, 

,'  .  -<z0-^z‘)lrz‘lm  (B-24) 

»;  -  <vzI*z?>WzSI.  (b-Z5) 

Rote  that  Z®,  Z®,  Z®,  V®,  and  v|  are  proportional  to  ZQ.  Hence,  they 
can  be  normalized  ty  Z^.  The  equivalent  network  ia  now  obtained. 


If  the  aperture  ia  symmetric  about  z*0  or  is  small,  then 


I  -  I 
1  NB 


NB  1 


(D-26) 

(D-27) 


By  (D-21)  -  (D-27),  the  network  is  reduced  to  the  symmetric  T 
network  shown  in  Fig.  D.7  where  the  impedances  and  sources  are 


z? 


z“ 


il'Wfa. 


2+f4*t 


NB 


2(,*WZ 0 


^l”*NB^2+V^NB^ 


^V^PV^nb 


f°  ■  (z  ♦Z®«-Z®)I  -Z®I 
2  '  0  1  r  NB  2  1 


(D-28) 

(D-29) 

(D-30) 

(D-30 


aacaaiSfiiM^^ 
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For  the  case  when  the  outward  traveling  TEM  currents  are 

symmetric  about  zm0  (I  *1  ,  I  *  I  _),  Z®  in  (D-29)  is  infinite  and 

1  V3  1  NB  2 

the  equivalent  network  in  Fig.  D.7  can  be  reduced  to  that  in  Fig.  D.8 
with  the  impedance  and  source  defined  as 

-2I„ Zrt 

Z®  •  - 1-2-  (D-32) 

V®  -  -I1(2Zq  ♦  Ze)  ( D— 33) 

They  agree  with  those  in  [4]  where  the  symmetric  case  was  considered. 
An  example  of  this  case  is  a  narrow  slot  whose  axis  is  perpendicular 
to  the  wire.  This  is  expected  because  the  magnetic  current  in  the 
slot  is  equivalent  to  a  superposition  of  magnetic  dipoles  which 
produce  symmetric  currents  on  the  wire. 


For  the  case  when  the  outward  traveling  TEM  currents  are 

A  A  Q 

antisymmetric  about  s*0  (!■-!.  I  «-I  ) ,  Z®  in  (D-2B)  vanishes  and 

I  NB  *  NB  • 

the  network  in  Fig.  D.7  is  reduced  to  that  in  Fie.  D.9,  where  the 
impedance  and  the  source  are 


<’-V2o 

A 

21. 


(D-34) 


Va  -  -(Z0+2Z®)l1 


(D-35) 
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An  example  of  this  case  is  a  saall  loop  aperture.  The  magnetic 
current  in  the  aperture  is  equivalent  to  an  electric  dipole  which 
points  normal  to  the  wire  and  produces  antisymmetric  currents  on  the 
wire. 


Fig. 


•  4* 


A  TEM  current  excitation  produces 

traveling  TRW  currants  f. e^kz«  and  I  e"1kzu  . 

1  ■“*  NB  -“a 


Pig.  D.5«  An  equivalent  T  network. 


outward 


Fig.  D.6.  Equivalent  circuit  for  the  calculation  of  V®  and  Ve. 


Fig.  D.7.  An  equivalent  symmetric  T  network  for  the  case  when  the 
aperture  is  symmetric  about  z-0  or  is  small. 


v# 


Fig.  D.8.  Equivalent  circuit  for  the  case  when  the  aperture  induces 
symmetric  outward  traveling  TEM  currents  on  the  wire. 


I  I 


z  «  0"  z  »  0* 

Fig.  D.9*  Equivalent  circuit  for  the  case  when  the  aperture  induces 
antisymmetric  outward  traveling  TEN  currents  on  the  wire. 


REFERENCES 


D.  Kajfez,  "Excitation  of  a  Terminated  TEM  Transmission  Line 
through  a  Small  Aperture,"  AFWL  Interaction  Note  215,  July 
t974. 

K.  S.  H.  Lee  and  F.  C.  Yang,  "A  Wire  Passing  by  a  Circular 
Aperture  in  an  Infinite  Ground  Plane,"  AFWL  Interaction  Note 
317,  Feb.  1977. 

S.  W.  Hsi  and  R.  F.  Harrington,  "Electromagnetic  Coupling  to 
an  Infinite  Wire  through  a  Small  Aperture  of  Arbitrary  Shape," 
Technical  Report  No.  19  on  Contract  No.  N00014-76-C-0225, 
Department  of  Electrical  and  Computer  Engineering,  Syracuse 
University,  Syracuse,  N.  Y.,  Feb.  1983. 

Yang  Naiheng  and  R.  F.  Harrington,  "Electromagnetic  Coupling 
to  an  Infinite  Wire  through  a  Slot  in  a  Conducting  Plane," 
Technical  Report  No.  15  on  Contract  No.  N00014-76-C-0225, 
Department  of  Electrical  and  Computer  Engineering,  Syracuse 
University,  Hatch  T982. 

C.  M.  Butler  and  K.  R.  Umasharikar,  "Electromagnetic  Excitation 
of  a  Wire  through  an  Aperture-Perforated  Conducting  Screen," 
IEEE  Trans.  Antennas  Propagat. ,  Vol.  AP-24,  No.  4,  pp. 
456-462,  July  1976. 

Yang  Naiheng  and  R.  F.  Harrington,  "Electromagnetic  Coupling 
to  and  from  a  Terminated  Wire  through  a  Rectangular  Slot  in  a 
Conducting  Screen,"  Report  No.  TR-82-7,  Department  of 
Electrical  and  Computer  Engineering,  Syracuse  University,  June 
4982. 

R.  F.  Harrington,  Time-Harmonic  Electromagnetic  Fields, 
McGraw-Hill  Book  Co.,  New  York,  1961. 

R.  F.  Harrington,  Field  Computation  by  Moment  Methods,  The 
Macmillan  Co. ,  New  York,  1968. 

A.  Sankar  and  T.  C.  Tong,  "Current  Computation  on  Complex 
Structures  by  Finite-Element  Method,"  Electronics  Letters, 
Vol.  11,  No.  20,  pp.  481-482,  Oct.  1975. 

J.  J.  H.  Wang,  "Numerical  Analysis  of  Three-Dimensional 
Arbitrarily-Shaped  Conducting  Scatterers  by  Trilateral  Surface 
Cell  Modelling,  Radio  Science,  Vol.  13,  No.  6,  pp.  947-952, 
Nov.  -Dec.  1978. 

J.  J.  H.  Wang  and  C.  Papanicolopulos,  "Surface  Patch  Modeling 
of  Scatterers  of  Arbitrary  Shapes,"  AP-S  Int.  Symp.  Digest, 
University  of  Washington,  Seattle,  pp,  159*16 2,  June  1979. 


1 15 


[l2]  S.  M.  Rao,  "Elect lomagnetic  Scattering  and  Radiation  of 
Arbitrarily-Shaped  Surfaces  by  Triangular  Patch  Modeling," 
Ph.  D.  Dissertation,  University  of  Mississippi,  August  1980. 

[ij]  C.  L.  I  and  R.  F.  Harrington,  "Electromagnetic  Transmission 

through  an  Aperture  of  Arbitrary  Shape  in  a  Conducting 
Screen,"  Technical  Report  No.  16  on  Contract 
No.  N00014-76-C-0225,  Department  of  Electrical  and  Computer 
Engineering,  Syracuse  University,  April  1982. 

[14]  A.  W.  Glisson,  "On  the  Development  of  Numerical  Techniques  for 
treating  Arbitrarily-shaped  Surfaces,  "  Ph.  D.  dissertation, 
University  of  Mississippi,  1978. 

[15]  R.  F.  Harrington,  "Matrix  Methods  for  Field  Problems," 
Proc.  IEEE,  Vol.  55,  No.  2,  pp.  136-149,  Feb.  1967. 

[46]  P.  C.  Hammer,  0.  J.  Marlowe  and  A.  H.  Stroud,  "Numerical 
Integration  over  Simplexes  and  Cones,"  Math.  Tables  Aids 
Comp.  ,  Vol.  10,  pp.  130-137,  1956. 

[17]  D.  R.  Wilton,  Private  Communication,  1983. 

[is]  J.  Van  Bladel,  Electromagnetic  Fields,  McGraw-Hill  Book  Co., 
New  York,  1964. 

[19]  H.  B.  Dwight,  Tables  of  Integrals  and  Other  Mathematical  Data, 
The  Macmillan  Co.,  New  York,  19^1. 

[20]  R.  E.  Collin,  Foundations  for  Microwave  Engineering, 
McGraw-Hill  Book  Co.,  New  York,  1^66. 


